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Abstract 

Color transparency (CT) in high momentum transfer (e, e'p) reactions is 
explored. The spin of the proton and photon are treated explicitly, hence the 
name "Vector CT". The Dirac distorted wave impulse approximation is used 
as a starting point; then CT effects are embedded. A hadronic basis is used to 
describe the struck proton as a wavepacket of physical baryon resonances. The 
effects of the wavepacket expansion on the normal component of the ejectile 
polarization, which vanishes in the limit of full CT, are investigated. This 
formalism is also applied to study CT effects in total cross sections, individual 
separated nuclear response functions, Fermi motion of the initial nucleon, 
non-zero size of the initial wavepacket and the effects of relativistic lower 
components. We show that including CT reduces the violations of current 
conservation (CC), a typical problem in calculations of this kind. The energy 
dependence of the normal polarization in (e, e'p) reactions is found to be 
slow. However, a measurement of the normal transverse response in a heavy 
nucleus, such as 208 Pb seems to afford the opportunity to see CT at quite 
low momentum transfers. The effects of Fermi motion are investigated, and 
choosing the momentum of the struck nucleon to be large leads to significant 
violations of CC. 

PACS numbers: 12.38Aw, 13.60-r, 13.85, 24.85. +p 
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I. INTRODUCTION 



Color Transparency (CT) is the postulated ]I],|2| absence of final- (or initial-) state inter- 
actions caused by the cancellation of color fields of a system of quarks and gluons with small 
spatial separation. For example, suppose an electron impinges on a nucleus knocking out 
a proton at high momentum transfer. The consequence of color transparency is that there 
is no exponential loss of flux as the ejected particle propagates through the nucleus. Thus, 
the usually "black" nucleus becomes transparent. We restrict our attention to processes for 
which the fundamental reaction is elastic, or at least a two-body reaction. This requires 
that the nuclear excitation energy be known well enough to ensure that no extra pions are 
created. This subject is under active experimental investigation 

The existence of color transparency depends on: (1) the formation of a small-sized wave- 
packet in a high momentum transfer reaction. (2) the interaction between such a small 
object and nucleons being suppressed (color neutrality or screening) and (3) the wave- 
packet escaping the nucleus while still small JI]^- That color neutrality (screening) causes 
the cross section of small-sized color singlet configurations with hadrons to be small was 



found in Refs. Ipf-Tlq], and is well-reviewed in Refs. JIB - So we take item (2) as given. 
The truth of item (1), for experimentally available energies, is an interesting issue. It is 
discussed in Refs. [|T§|-f2~T|] , and is not probed in depth here. 

It is also true that, for present experiments PHlOU, the small object does expand as it 
moves through the nucleus. Thus the final-state interactions are suppressed but not zero. 
The importance of this expansion was found by Farrar et al. |2~2[, and by Jennings and 



Miller p^-g6j|. See also Refs. |3,g8 



Up until now, all calculations of CT have totally neglected the effects of spin. In this 
paper, the full Lorentz structure of the matrix element for the electroproduction of nucleons 
from nuclei is considered. We call this work Vector Color Transparency, where the vector 
moniker derives from the vector nature of the photon exchanged in the hard collision. 

What do we hope to gain by introducing this complication? First, experience has shown 
that the introduction of spin into many physical problems led to totally new phenomena. 
Whether or not this happens in discussions of CT is the subject of this paper. Second, by 
invoking time-reversal and parity invariance of the electromagnetic and strong interactions, 
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the normal component of the polarization of the ejected proton [ in (e, e'p) experiments] must 



vanish in the limit where final-state interactions are absent [25]. This is the limit of full CT. 
Therefore, the spin dependent observables provide a very sensitive measure of the effects of 
CT. Thirdly, we will soon see that the spin dependent observables are defined in terms of 
ratios of sums of response functions. Thus, experimental measurement of such observables 
are less susceptible to systematic errors than, for instance, the unpolarized observables. 
Lastly, an experiment has been proposed and approved, to be run at CEBAF ||, which will 
measure the energy dependence of the unpolarized cross section and the normal component 
of the ejected proton's polarization. In that case, it is desirable to have ready theoretical 
predictions, crude though they might be, to be confronted by the experimental results. 

In a previous paper |3(| a multiple-scattering series for the color transparency 
wavepacket-nucleon interaction, within the framework of Glauber (eikonal) theory, was de- 
veloped. The approach here is the same, except that four-component Dirac spinors are used 
for both scattered wavefunctions and bound states. As in Ref . , we embed CT operators 
and states into the usual DWBA treatment. 

The organization of this paper is as follows. In Section II we display formulae for the 
(e, e'p) cross section and polarization in terms of the nuclear response functions, the hadronic 
tensor, and the nuclear current matrix element (NCME). In Section III the scattering state 
wavefunction is derived within the context of the Dirac impulse approximation. Most of the 
material in Sections II and III is well-known, but we present these sections as necessary pre- 
liminaries to including the effects of CT. This is done in Section |VL by including the effects 



of the composite nature of the nucleon. We develop approximation schemes, as in Ref. [30 



to evaluate the CT wavefunctions, in Section M. Section VI contains a demonstration that 
there is at least one reliable approximation scheme for each range of energies. In this section, 
we also present numerical results for total cross section ratios, integrated longitudinal and 
transverse nuclear responses, current conservation violations, differential unpolarized cross 
sections, differential normal polarizations, individual nuclear response functions, effects of 
non-zero wavepacket size, Fermi motion, and the effects of the purely relativistic lower com- 
ponents of the bound and scattered state wavefunctions. In the final section, Section VII, 
we summarize and make some concluding comments. 
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The publication of the results of the SLAC experiment || is imminent but not published 
now. It is important to use the experimental acceptance in computing observables. Thus 
we leave a detailed assessment of that experiment to a future publication. We also make no 
attempt to completely review the CT literature. 



II. THE NUCLEAR CURRENT MATRIX ELEMENT AND CROSS SECTION 



We now describe the spin-dependent formalism. In particular, consider the case of bom- 
barding a spin-0 nucleus with polarized electrons and detecting polarized protons which are 
knocked out. This process is usually written as (e,e'p). The incoming electron has helicity 
h and the outgoing nucleon is polarized, in its rest frame, along sr. A schematic drawing of 
the (e, e'p) reaction is shown in Figure [I]. The following discussion reproduces the essential 
points of Ref. BT] 



The kinematical situation is shown in Figure |2|. We take the virtual photon to lie along 
the Z direction, and the electron scattering plane to be the YZ plane. The scattering 
amplitude for Vector CT is defined as 

M a = j ll J£, (2.1) 

where is the matrix element of the electron current and J£ is the matrix element of the 
nuclear current (NCME). In the one-photon exchange approximation and using the nuclear 
shell model, the cross section is 

d s a _ M N \p\ f da \ 



Matt 



(2tt) 3 \dn k - j 

x £ [dE p \M a \ 2 5(E p -q°-M N + e a ), (2.2) 

where k^(k'^) is the four-momentum of the incoming (outgoing) electron, = — k" 1 . 
The binding energy of the nucleon in shell a is e a . This is small compared with other 
energy scales in the problem and is neglected. The sum is over all occupied shells. The 
four- momentum of the outgoing nucleon is p M and its energy is denoted by p° = E p ; M N is 
the nucleon mass. The solid angle dQ p = sin £ d( d/3. The nucleus is assumed to recoil with 



negligible energy. The Mott cross section is defined in Ref. [31 
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The square of the matrix element M a can be written in terms of electron and nuclear 
tensors: 



\M a \ 2 = r)^Wr, (2.3) 

The electron tensor r]^ v = | (k^k '„ + k v k' ^ — g^k ■ k' + iht^ v x K k' x k K ^ where is the met- 
ric tensor with g 00 = —g\\ = — #22 = ~933 = 1 ; an d is the completely antisymmetric 
fourth rank tensor. The electron helicity is given by h. Here and throughout we use the 
notation and conventions of Bjorken and Drell [|32] . 



We turn now to the nuclear tensor. This depends on the direction of the (rest-frame) spin 
of the ejected proton. Thus, W£ v = W£ u (sr) where we use the notation that a describes a 
unit vector in the direction of a. The upside-down "hat" is used to distinguish unit vectors 
from operators. The nuclear tensor is 

WHS*) = J£*(q)JM, (2.4) 

where J^(q) are the matrix elements of the nuclear electromagnetic current operator. All of 
the unknown physics is lumped into the nuclear current. Therefore, we define the quantity, 

(S fl ) =Y,i dE P W « U is*) S ( e p -q°~M N + e a ) , (2.5) 

where the sum is over all occupied single-particle shells and the integral and delta function 
serve to conserve energy in this full nuclear tensor. The nuclear tensor W£ u depends on the 
ejected proton energy, E p , because the scattered proton is ultimately detected, on the mass 
shell, moving with momentum p. 

Thus, the crucial quantity is the nuclear current matrix element (NCME). We now 
present a symbolic form for this current. The NCME is given by 

JM = { -\%M(q)\a). (2.6) 

where the initial state of this knockout process is labelled by the shell model state of a bound 
proton, a. In the final state, a proton moves with a momentum p with some rest frame spin 
Sr. The residual nucleus is an a-hole state, which is not detected and is therefore summed 
over. The initial and final states are connected by a vector operator denoted as T^(q) which 
describes the absorption of a virtual photon on a proton in the nucleus. The overlap of the 



initial and final nuclear states is imagined to be a single particle state of a nucleon bound in 
shell model state a. Incoming boundary conditions for the proton scattered wave are used, 
but we can also use outgoing boundary conditions: 

= m\q)\%,-s R ) {+) ■ (2.7) 
A more complete derivation is described elsewhere BT] . 



We contract the indices of the electron and nuclear tensors to obtain the differential cross 
section. The use of current conservation: q^J^ = q°J° — q 3 J 3 = 0, where q = q 3 Z allows us 
to eliminate the longitudinal components of the nuclear current in favor of the better known 



2 

charge density. In particular, W 33 = ^W 00 and W 3k = ^W ok for k ^ 3. We note here an 
important point concerning current conservation and gauge invariance. The NCME is an 
integral over a bound state wavefunction, a scattered state wavefunction and an operator 
which connects the two. If these two wavefunctions are derived from the same theory, as 
they are in nature, then as members of a complete set they are orthogonal and current is 
conserved. If the bound and scattered states are modelled with different Hamiltonia, as 
usually practiced, there is no guarantee that the wavefunctions are orthogonal and that 
current is conserved. This problem is discussed in Section |V] where we show that current is 
conserved in our calculation at the < 10% level or better. 

The above discussion leads to the differential cross section for the scattering of polarized 
electrons (helicity h) off of a nucleus, ejecting a proton with rest- frame spin s# pi] : 



, x / d 3 a \ 
*h{sR) = J ^ ^ (2.8) 



^dek'dVtk'dVL p j h ^ 
Mjv|p| / da 
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{V L (R L + R n L S n ) + V T (R T + R n T S n ) 



2(2tt) 3 \dtt k , 

+V tt [(Rtt + RttSu) cos 2(3 + (R l TT S t + R l TT S t ) sin 2(3} 
+V LT [{R LT + R n LT S n ) sin (3 + (R LT S t + R l LT S t ) cos/5] 
+hVl T [(R' LT + R'2 T S n ) cos (3 + (R!l T S t + R a LT S^) sin (3} 
+hV^ T (^R'r^ T St + R'tt<~>i)}i (2-9) 

where kinematical factors V are defined by Vl = {fa) > = tan 2 \ — ^ TT = ~ 2jqP' 
V LT = fa (tan 2 6 f~fa) , V' LT = fa tan f , V± T = tan f (tan 2 f - ^ 



is the electron scattering angle, q 2 = q^, and |q| is the magnitude of the photon three- 
momentum. The direction of the ejected proton's spin is given, in its rest frame, by Sr. The 
spin projections which appear in the above cross section are identified to be: S n = n • Sr, 
Si = 1 • Sr, S t = t • Sr, where the unit vectors are n = (3 = (q x p) /|q x p|, 1 = p 
and t = ( = nx 1 . The quantities R are called the response functions. The subscripts 
L and T are the contributions from longitudinal and transverse photons (both helicities), 
respectively. The subscripts TT and LT refer to the transverse- transverse and longitudinal- 
transverse response of the system. These arise from the interference between the different 
polarizations of the photon. The subscripts TT' and LT' also arise from the interference of 
the longitudinal and transverse photons. They are "primed" because they are formed from 
an antisymmetric combination of components of the nuclear tensor and thus only accessible 
with a polarized electron beam. The response functions themselves are defined in terms of 
the hadronic tensor. Explicitly, 

l -{R L + R n L S n ) = w"\s R ), (2.10) 

l - (Rt + R T S n ) = W n (s R ) + W 22 (s R ), (2.11) 

^[(R TT + R TT S n ) cos 2/3 + (R^St + R l TT S t ) sin 2/3] = W 2 \s R ) - W U (s R ), (2.12) 

^[(R LT + R n LT S n ) sin/3 + (R^St + R 1 ^) cos/3] = W° 2 (s R ) + W 2 °(s R ), (2.13) 

1 -[{R' LT + R'2 T S n ) cos/3 + (R'l T S t + R' l LT S) sin/3] = % [W 10 (s R ) - W°\s R )} , (2.14) 

\{{R^rS t + R* T Sty = % [W 12 (s R ) - W 2 \s R )] . (2.15) 

Since the response functions, as defined above, are independent of /3, certain restrictions 
can be imposed using simple kinematical arguments. In particular, in "parallel kinematics" 
when p = q, all of the response functions must vanish except for Rl, Rt and R'^ T , which are 
unconstrained, and the normal and transverse parts of Rlt and R' LT which are constrained 



to be opposites |3T| , |33|1 . We point out that the constraints of parallel kinematics have been 



violated in the work of Allder [34]. Refs. [31 1 and [33] show how these response functions 



can be extracted from the cross section. 

Now, we leave the general theory of the (e, e'p) reaction and focus on the upcoming 
CEBAF experiment ||. Saha and collaborators plan on measuring the unpolarized cross 
section and the normal component (to the photonuclear scattering plane) of the proton 



polarization. In terms of the above notation, the unpolarized differential cross section is [35 



+V TT R TT cos 2f3 + V lt Rlt sin (3} } (2.16) 

where the subscript on a signifies that electron helicity is averaged and the argument of 
implies that the direction of the ejected proton spin is averaged over. The normal component 
of the polarization vector is then defined as 

o _ 1 M N \p\ / da \ 

ct (0) {2ny \dil k 'J Mott 

+V TT R^ T cos2p + V LT Rl T smp}. (2.17) 

Since the normal component of the polarization, P°, is dependent only on normal response 
functions which vanish in the high energy of the CT limit, P° also vanishes in that limit. 
Parity and time reversal invariance imply that non-zero values of P° arise only from the 
final-state spin-orbit interaction. If CT is obtained all final-state interactions vanish, so 
that P° vanishes too. 

It is also useful to consider the "total" unpolarized cross section, obtained by integrating 
over the solid angle of the outgoing proton. If we integrate over (3 in Eq. ( |2.16|) , all of the 
terms which have explicit (3 dependence vanish. It then remains to integrate over sin(d(. 
However, at this stage, we introduce the convenient approximation that dfl p ~ d 2 \p— q|/|q| 2 . 
Since the response functions are functions only of q^ = p — q, (q± = |q_i_|) we define the 
total unpolarized cross section as 

_ d 2 a ^ ^ 

da^idVly 
Mat|p| I da \ 
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(2ny |q| 2 \dn k , 

x J dq ± q ± {V L R L (q ± ) + V T R T (q ± )} . (2.19) 



S 



The above equations completely specify the relationship between the response functions 
of the (e, e'p) reaction and the experimental observables which will be measured in the near 
future. The crucial quantity is the electromagnetic nuclear current matrix element (NCME) 
of Eq. Q. 

III. DWBA FORMALISM 
A. Bound State 



We now discuss the model bound state used here. In Ref. 0], the ejected proton 
was taken to lie initially in a non-relativistic one-particle harmonic oscillator shell-model 
state. Since the present formalism is spin dependent and relativistic, it is necessary to look 
for a relativistic single-particle shell-model state which is a four- component Dirac spinor. 
There are several such models |BB|,|57|1. Here we use the finite nuclei model of Horowitz and 



Serot Iftql - Our opinion is that QHD is not a fundamental field theory. However, for our 
purposes it is not necessary to consider QHD to be anything but a model, which accounts 
for nuclear phenomenology. In particular, QHD yields realistic nuclear densities and four- 
component bound state wavefunctions. 

It is useful to display the solution to the Dirac equation of QHD in terms of radial and 
angular solutions. The set of quantum numbers a = {n, j,l,m,t} where n is the principal 
quantum number, j is the total angular momentum, I is the orbital angular momentum, m 
is the Z component of the total angular momentum and t is the third component of the 
isospin. Then, the solution is 



(R|a) = $ Q (R 



1 
R 



iG n ji(R) 
+F n ji(R)er ■ R 



yjlm{R)Vt, (3.1) 



where G and F are functions of the radial variable R — |R| only and rj t is a two-component 
isospinor, with r\ t = | for protons and r\ t = — | for neutrons. The quantities [Vj/m are the 



standard spin spherical harmonics. We use a computer program by Horowitz |39| to generate 
these radial wavefunctions. Note also that $ Q (R) are normalized to unity. 
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B. Current Operator 



The vector current operator represents the absorption of a virtual photon by a nucleon 
bound in a nucleus. This nucleon is somewhat off the energy shell. In general the current 
operator can be expressed in terms of scalar functions of the four-momenta multiplied by 
any of twelve four- vectors [|3B| ! Since the bound nucleon is only slightly off shell (the binding 
energy is small), and since any other choice for this operator is only a guess, the current 
operator is chosen to be the free nucleon current operator. Thus, we define the configuration- 
space matrix elements, which are matrices in the Dirac space, as, 

(R|T£(g)|R'> = e iqR (5 3 (R-R')7°r^(g), (3.2) 

where 

r"(g) = Fi(?V + ^^(? 2 K^, (3-3) 

and Fi(q 2 ) and F2(q 2 ) are the usual Dirac and Pauli form factors of the nucleon. We use a 
simple dipole parameterization such that 

^M'-imM" (3 ' 4) 

G M {f) = (1 + k)G £ (<j 2 ), (3.5) 

where k — 1.79 is the anomalous magnetic moment of the proton. The Dirac and Pauli form 
factors are given by 

Fi{q2) = t G M (£-^G i <g) t 



q 2 - AM 



N 



2 

m 2 ) = q2 T 4 N M 2 N [ G ^ 2 ) - g m^ 2 )\ ■ (3.7) 

This assumption of the current operator neglects the possibility that the nucleon properties 
are much modified in medium. 



4M 



C. Optical Potential and Distorted Wave 



We now describe the wavefunction of the scattered proton. Here the bound state proton 
wavefunction is a relativistic four-component object, and the current operator is a 4 x 4 
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matrix. Thus the proton wave must also be a Dirac spinor, satisfying a one-particle Dirac 
equation. 

Treating the Dirac equation as a one-particle equation has many serious diseases. Indeed, 
problems of interpretation with the negative energy states led to the development of quantum 
field theory. However, the difficulties can be identified (see e.g. Section 2 of Ref. JKJ.) Here 
we are interested in the scattering of energy eigenstates of an interacting Dirac Hamiltonian. 
Such eigenstates propagate along with no transitions to the negative energy states. In these 
situations, there is no problem in the interpretation of the Dirac equation as the governing 
equation of a single particle relativistic quantum mechanics. 

In the early eighties, a phenomenology based on using the Dirac equation to describe 
the scattering of medium energy protons from nuclei was developed. Both partial wave 
analyses |4l| , |42[| and eikonal reductions of the Dirac equation []43|-|46f have been performed. 



The successes of these models, most notably the spin observables such as analyzing powers 
and spin rotation functions, as well as differential cross sections, are still quite impressive. 
From the standpoint of CT, where we are interested in very high energy scattering, the 
eikonal form is appropriate. 

Let us now discuss the potentials that enter into the Dirac equation. In the impulse 
approximation these are complex Lorentz scalars and time- components of 4-vectors. We 
label these V s and V v , respectively. In the relativistic impulse approximation, the scalar 
and vector optical potentials, V s and V v , are proportional to forward Dirac scalar and Dirac 
vector scattering amplitudes (F®, F®) and to the scalar and vector nuclear densities [f|7]-|53[ 
(p s , Pv)i such that 

U opt (R) = V S (R) + 7°K(R), (3.8) 

where 

V 3 (R) = rF°p a (R) and V V (R) = rF°p v (TL), (3.9) 

and r = —ATripi^/Mx is a kinematical factor needed to relate the invariant Feynman am- 
plitude to the usual description of the optical potential in the impulse approximation. The 



scalar and vector nuclear densities are taken from the QHD model in the previous section [38 . 
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Ps(R) = E 



2jg + 1 

AttR 2 



(\G a (R)\ 2 - \F a (R)\ [ 



= £ ) (\G a (R)\ 2 + \F a (R)\ 2 



(3.10) 
(3.11) 



The potential strengths, _F S ° and F® are taken directly from NN scattering data [pO|,^-|6T 
see Appendix A. Since the strengths are a sensitive function of energy and NN scattering 
data, we consider our calculations reliable only at the energies where sufficient data exists 
to determine the strengths F® and F®. 

The Dirac equation for the distorted wave is given by 



-icx-V + P(M N + V s (R)) 

+V v (R)W+l(R) 



(+) 



(3.12) 



where E = Jp 2 + Mjy. This equation is solved by separation into two coupled first-order 



differential equations. The eikonal form of the solution to this equation is well known [fl3-46 
and is given by 



1 

icr-V 



x exp 



E+M N +V 3 (R)-V V (R) 
Z 

dZ'Q(B, Z') 



X-s 



(3.13) 



where R = B + ZZ, 

Q(B,Z') = ^-{U C (B,Z') 

lip l 



+U so (B,Z') [a-Bxp-tpZ']}, 



(3.14) 



R ■ p = pZ, and B ■ p = 0. Also note that we have used outgoing boundary conditions and 
changed the sign of the rest-frame spin projection, as in Eq. ( |2.7| ). The normalization is 
J\f = J(E + Mn)/(2Mn). The Dirac scalar and vector potentials have been eliminated in 



favor of new central (U c ) and spin-orbit (U so ) potentials: 

U C (R) = 2EV V (R) + 2M N V S (R) 

+V S 2 (R)~V V 2 (R), 



(3.15) 
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E + M N + V S {R)-V V (R) 

x R~M {Vv{R) ~ Vsm - (3 - 16) 

Note that the potentials U c and U so have units of {Energy) 2 ] we use the symbol U (and not 
V) to indicate this. The expression in Eq. fl3.13| ) can be simplified; the detailed form of the 
wavef unctions is relegated to Appendix B. 

Finally, we note one other complication. The eikonal wavefunctions presented in Ap- 
pendix [FJ have been derived in the limit where the particle travels along the Z axis (direc- 
tion of the virtual photon). At very high energies, the difference between p and Z is small. 
Naively, this finite scattering angle should make little difference. However, if one assumes 
that the scattering angle is exactly zero, then the response function Rtt is also exactly 
zero. Since, in the plane-wave limit, we know this is not true, account must be taken of this 
scattering angle. We have derived the cross section in Eq. (|2.9| ) with the coordinates that 
the final proton momentum makes an angle ( with the Z-axis (the direction of q) and an 
angle (3 with normal to the electron scattering plane (the XZ plane). Call this the g-basis. 
The derivation of the eikonal wavefunction described above has the proton momentum p 
along Z; call this the p-basis. To move our eikonal wavefunctions from the p-basis to the 
g-basis, we perform a passive rotation of the coordinates to the desired frame. 



IV. INCLUSION OF CT 



To include the effects of color transparency, we proceed as in Ref. JHJ and include the 
baryonic internal degrees of freedom. Thus, the states in the full Hilbert space now carry 
two labels: one for the "external" nuclear space and one for the internal quark space. The 
quark space operators are identified wearing "hats" and we write CT quantities using script 
lettering, where possible. Then, the CT scattering matrix element is given by 

M a =j,JZ (4.1) 



where the electron current is j^. The triple differential cross section is given by Eq. ( j2.2|) , 
except for the replacement of M a by Ai a . Thus, in order to completely specify the method 
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for the inclusion of CT effects, it remains to construct the CT generalized NCME. We follow 
the method in Ref. |5D[1 and write 



J£(q) = (N,a\fg(q)\* N>p ^ R )M, (4.2) 

Here, |^at, p ,-s h )^ + ' ) is a vector in both internal and external spaces. The subscript on 
\l/ is to indicate that ultimately a nucleon N moving with momentum p and rest-frame 
spin projection — sr is detected. The overlap of this state with the nuclear position state 
(R| gives a vector in the internal space only, which we identify with bold face type: 
(R|\I / Ar i p i _s ii )^ + ' ) = ^tv~p _g B (R)- The asymptotic boundary condition is implemented by 
imagining that Vf^p,-^ is some internal space vector, whose equation of motion we shall 
derive below, projected onto the nucleon state \N). 

In our approach the relativistic bound state is not changed by including CT effects. See, 
however, Ref. |62[ . The rest of this Section describes the other elements in Eq. ( |4.2j ). 

A. Current Operator 

Upon absorption of a high-energy photon, the proton is converted into a coherent super- 
position of baryon states (a wavepacket), with components labelled by a discrete quantum 
number, m. Then, we identify the relevant inelastic transition Dirac and Pauli form factors 
by taking the matrix elements of the current operator to be 

(iV,R|f^(g)|m,R') =e^ R 5 3 (R-R')7°r^ m (g), (4.3) 

where 

B. Optical Potential 

We now consider the optical potentials V S (R) and V V (R) to be operators in the internal 
quark space. These optical potentials are essentially products of the baryon-nucleon scat- 
tering amplitude operator f{b 2 ) with p s or p v . We choose the scalar and vector potentials 
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to have the same operator structure, to be proportional to the same function f(b 2 ). At 
present, there is no detailed knowledge about the precise form for f(b 2 ). However, some 
general properties are known. For small wavepackets with b 2 <C b 2 H with b 2 H = (N\b 2 \N), 
f should vanish. Interactions do occur for larger wavepackets. For non-zero but still small 
sized wavepackets, the interaction goes like 

b 2 



Jim f(b 2 ) -> w . (4.5) 
The operator / is normalized such that 



f{b 2 = 0) = 0, (4.6) 
(iV|/(6 2 )|A0 = l. (4.7) 

We use a general function f{b 2 ), subject only to the constraints of Eqs. ( |4.5|) - (|4~7|) . The 
matrix elements of the scalar and vector operators can then be written as 

<ro,R|V 8(t0 |m',R') = 5 3 (R-R')V s{v) (R)(m\f(b 2 )\m'}. (4.8) 

To compactify the notation, we define the quantity V s ( v )(R)f(b 2 ) = V s r v )(H). 

C. Distorted Wave 

In this section we obtain the wave equation for the propagation of the wavepacket formed 
in the hard collision through the nucleus. Start by considering the time-independent Dirac 
equation. The equivalent Dirac Hamiltonian should now be considered an operator in the 
internal quark space. As such, we award H a "hat". All of the tools have already been 
developed and we immediately write 

W+l iSh (R) = [-ia ■ V + /3 (M + V s {R) 

(+) 



+V V (R) ^^(R) 

N,p,s R \~ D ' 



E9™(R), (4.9) 



where E is the energy and we have already taken the matrix element in the external con- 
figuration space. That is, ^^ ^(R) = (R|^Ar )Pi s fl ) ( - + ' ) . In this notation, M 2 is the baryon 
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mass operator squared. That is, M 2 \m) = M^Jm); the nucleoli is the ground state with 
m = N having eigenvalue Mj^. The combination 



P 



2 = E 2 - M 2 



(4.10) 



is also a quark space operator which can be interpreted as the baryon momentum operator 
squared. This operator accounts for the different kinematics with which the wavepacket 
components propagate through the nucleus. Here p is also a three- vector in the Z direction; 
p = pZ. 

The Dirac equation can again be solved by eliminating the lower components. Then 
we see that the CT scalar and vector Dirac potentials (which are operators in the internal 
space) can be eliminated in favor of combinations given by 



U C (R) = 2EV V (R) + {M, V s (R)} + [M, V v (R) 



+V 2 (R)-V 2 (R), 



Uso(R) = ^ 



^(y v (R)-V s (R) 



(e + M + V s (R)-V v (R)) \ 



(4.11) 



(4.12) 



where the square (curly) brackets are the (anti-) commutator symbols. We also define 
"path-evolved" versions of the operators of Eqs. fl4.11|) and ( (4.12| ) as 



U c{so) {K)=e-^ z Jj c{so) {R)e 



ipZ 



(4.13) 



Eikonalizing in the usual way, we see that the solution to the resulting eikonal first-order 
equation is a path-ordered exponential. The notation is simplified by defining an operator 
Q such that 



1 
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U C (B,Z') 



+W so (B, Z') {a ■ B x p - ipZ') 



(4.14) 



Note that 
and 



f2(R), f2(R') 7^ since p and b do not commute. This is because p contains M 2 , 
M 2 ,f(b 2 )}^0. 

We want to identify the equivalent CT wavefunctions which are obtained when the NCME 
is written as 
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JM = J d 3 R<S> a (R) e-^T^q) tf$ iP ,_ 8jI (R). (4.15) 

where $«(R) = (R)7° and F = 7 [r^] f 7 . The bound state, $ a (R), is given in Eq. (PD 

Quite generally, then, the CT 



and the nucleon current operator is defined in Eq. (B.2| 
wavef unction is given by 



*S,p,-^( R ) = ^E[H?) 



{ml 



( 



V 



1 

icr-V 




dZ'Q(B,Z') 



\N) X - 



S R 1 



(4.16) 



S+M+V S (R)-K(R) 

where we have taken the momentum of the outgoing wavepacket to lie along the Z direction 



in Eq. ( |4.16| ) but not in the calculations. Again, M = y (E + Mjy) / (2M^). No sum over the 
index /1 is implied in the above equation. Above, V is the path-ordering symbol. The effects 
of CT effects can be seen by comparing this expression with Eq. ( |3.13| ). We have projected 
the column vector St on to the nucleon state to be detected, and could also project on to a 
nucleonic isobar to obtain the isobar production amplitude. The inelastic nucleon current 
operator, r^ m , is defined in Eq. ( |4.3| ). 

The expression for the CT wavefunction, shown in Eq. ( |4.16| ), is the central formal result 
of this paper. Note that we have changed the direction of the rest-frame spin projection, in 
accordance with Eq. 



1. Approximations: Vector OBO, LEE and EA 



At this stage, we may proceed to make the expansions and approximations, as in Ref . |30| . 
which allow a numerical evaluation of the CT wavef unctions, Eq. ( [1.16j ), the NCME, 
Eq. (|4.15| ) and ultimately of the (e, e'p) cross section. The definitions of these approxi- 



mations are only slight generalizations of those in Ref. p0| |. The "Order by Order" (OBO) 



scheme is defined by simply expanding the path-ordered exponential in Eq. ([4.16|) order by 
order in its Taylor series and evaluating the resulting matrix elements between the PLC 
and the nucleon. The Low-energy Expansion (LEE) scheme is defined by factoring out the 
expectation value of this operator Q in the nucleon by letting f2(R) = (iV|f2(R)|iV)-|-Af2(R) 
and treating AQ as a perturbation. One must be a little careful, because of the presence of 
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the noncommuting Pauli matrices. We have also calculated the "Exponential Approxima- 
tion" (EA) defined by exponentiating the first-order result of the OBO (this is equivalent 
to neglecting the path ordering). This is similar to letting 1 — x ~ e~ x . Our previous 
experience leads us to believe that the EA and the LEE will be accurate at the low and 
medium energies while the EA and the OBO should be reliable at the highest energies. This 
is verified below. 



V. MODEL EVALUATIONS AND APPLICATIONS 

A. Wavepacket-Nucleon Interaction and Quark Space 

To proceed further we need specific forms for / and the states m. We choose the inter- 
action to be 

f(P)=P/b 2 H , (5.1) 



as in Eq. (4.5). Further, we choose the internal baryonic states to be described by a two- 
dimensional transverse harmonic oscillator. For interactions of the form f(b 2 ), the use 
of three dimensional oscillators leads to the same results as those of the two-dimensional 
case. We calculate transparencies for bn = 1 fm and for two different values of the oscillator 
spacing: Mf — M% = 1.19 GeV 2 and Mf — Mjj = 2.36 GeV 2 . These simple choices reasonably 
represent the present (lack of detailed) knowledge of f(b) and the baryon wavef unctions. 
With the internal model space and interaction operator specified, we note the useful formula: 

P 
9 



( 2m ll2~l 2n ) = ( 2m + X ) <W -m8m,n+l 



-(m + 1) 5 m+1>n . (5.2) 

Before we describe our evaluations, it is necessary to show the parameters we use for 
the strengths of the optical potentials. These are determined in Appendix |A| and displayed 
in Table |I[ It is useful, when examining the table, to recall Eqs. ( |3.8| ) and ( |3.9| ). Also, 
recall Eqs. ( |3.15 ), ( |3.16| ), (4.11) and ( |4.12| ) to see that at high energies the central potential 



dominates and that, within the central potential, the vector potential is the important 
quantity. Thus, we can immediately see that at high energies the optical potentials we 
display in Table | are absorptive as the usual optical potentials. 



The last column of Table | shows the potential strengths in nuclear matter in more 
conventional units, where we have taken p = 0.166 fm~ 3 . The optical potential strengths 
are energy dependent, but this energy dependence is unknown for higher energies where 
there is no data. Therefore, we have no way to really know what the strengths should be 
at higher energies. The strengths are totally determined by NN scattering data, but they 
are functions of more than just the invariant cross sections, see Appendix [A[ Thus, we only 
consider our calculation to be reliable at energies where enough data exists to at least allow 
an approximate determination of the strengths. 

As described in Appendix |A], the strengths displayed in Table | can be arranged to give 
the proton-nucleon total cross section and the ratio of the real to imaginary parts of the 
forward scattering amplitude, a/. These numbers are summarized in Table [II]. 

Defining some more notation helps us proceed. First, we introduce the shorthand nota- 
tion U^'fi such that 

(2m\U c (R)\2j) = U {2m > 2l) (R) (5.3) 
(2m\U so {R)\2j) = U?™*\R). (5.4) 

In the above equation, it is understood that JJ^ ±2n ' 2 ^' vanishes for n > 3 while JJ < ^ ±2n ' 2 '^ = 
for n > 2. 

More explicitly, we can evaluate the matrix elements just defined in terms of the scalar 
and vector optical potentials, using Eq. (|5.2| ). Evaluation of the spin-orbit potentials require 
a little more discussion. In particular, the inverse operator appearing in Eq. (|4.12|) is prob- 



lematical because 



M, V s{v) 



7^ 0. However, neglecting the operators K(R) and V V (EL) in 
this inverse operator is of the same order of approximation as those already made in deriving 
the eikonal approximation so that it is safe to replace 



-i 



-l 



E + M + V S (R) -K(R)J ~{E + M\ , (5.5) 

also in Eq. ( |4 . 1 6| ) to be consistent. In that case, we can obtain explicit representations for 
the spin-orbit CT potentials. As an example, 

«>^^^ 
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The first few central CT potentials have the explicit representations 

U^' 0) (R) = 2EV V (R) + 2M N V S (R) 

+2V S 2 (R)-2V 2 (R), (5.8) 
U?>°\R) = (M N -M 2 - 2E) V V (R) 

- (M N + M 2 ) V S (R) - W 2 (R) + W 2 (R), (5.9) 
U^\R) = 2V 2 (R)-2V, 2 (R). (5.10) 

where we have only displayed the nonzero CT potentials which are connected to the ground 
state (nucleon). Since only the spin-orbit terms enter with the spin operator, it is convenient 
to define 

U {2 T 2j \K) = Uf m ' 23 \R) - ipZ U {2m ' 2 ' J \R). (5.11) 

It is also useful to define, here, a commonly appearing combination of the above CT 
potentials: 

oo 

<£)(B, Z, Z>) = £ U%*»(B, Z)e^-^ z - z '\ (5.12) 

m=0 

In performing the calculations, one also encounters terms like -fj- and which appear 

inside the summation sign above. In the spirit of the high-energy nature and approximations 
of this paper, in the calculations we assume all of these quantities to be equal. Thus, the 



expressions appearing in Eq. (|5.6|) and Eq. (|5.7|) are opposites in this approximation. The 



differences introduced by these approximations are smaller than the other uncertainties in 
the calculation. 



B. Dirac Effects 



We point out some interesting results obtained with this Dirac formalism. This is in- 
cluded in anticipation of the numerical results discussed below. The following effects are 
general features of a relativistic description of the wavepacket. 

1) The presence of the commutator term in Eq. ( [4.1 1|) causes a new non-Hermiticity in 
U C (R). Thus, care must be taken in order to obtain the correct sign of the commutator term. 
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Because of the quark operator structure present in U c , and in U so , it is a little bit tricky 
to switch between outgoing wave and incoming wave boundary conditions. In the DWBA, 
all we do is take the complex conjugate of the potentials. Here, that will not suffice. One 
must take the complex conjugate of the strengths, while leaving the CT part alone, except 
for making the change in this commutator term. Because this switch is complicated and 
relatively subtle, we only consider outgoing wave boundary conditions. We have checked, 
however, that both descriptions give the same result. Hereafter, this equality is to be 
understood and we drop the superscripts indicating the boundary condition. 

2) The presence of the quadratic terms in the central potential causes an inequality 
between quantities UfP\R), as defined in Eq. ( |4.11| ), and U C (R) as defined in Eq. (P.15Q . 
This effect is an artifact of the way we solve the Dirac eikonal formalism. 

3) The commutator term in U c alters the approach to full transparency. This is seen in 
the following way. In the limit of zero size wavepacket, the inelastic and elastic form factors 
are equal. Then, we can neglect the subscripts on iq, F?, and r M . In that case, the initial 
PLC is described simply by \PLC) = f£(q)\N) = T^q) Em\ m )- Then > 

(PLC\U C (R)\N) = r»(q) ^2(m\U c (R)\N), (5.13) 

m 

= r"(?)AM (V S (R) + V V {R)) , (5.14) 

where AM = Mjy — M 2 . The reader is urged not to panic, however, despite the appearance 
that in the closure limit CT is not obtained. Upon looking at the definition of the CT 
potentials which actually enter, Eq. ( 4.13 ), we see that U c always enters with a factor of 



p in the denominator. Thus, in the high-energy limit (closure limit) p ~ p — > 00, the 
exponential factors do not matter, and the i factor will suppress Eq. ( |5.14j ). This is not 



unlike the situation in Ref. ||30|| . In that case, the approach to transparency was described 
as 1 — e* ( - P2_p - )A ~ (p2 — p)A for large energies. Since p ~ p + (Mjy — M 2 )/(2p), we see that 
Pi—p~ (Mjv + M 2 )(Mjv — M 2 )/p. So, we see that CT is still obtained at very high energies, 
only its approach is affected. Furthermore, the terms which are left over, in the CT limit, 
are of the same order as terms like p/p 2 j which we have neglected already. 

4) New color transparency derivative (CTD) effects enter. To illustrate this effect, con- 
sider the one-dimensional functions 
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F(Z)= f Z dZ'f(Z'), (5.15) 
F CT (Z) = f Z dZ'f(Z') \l - e *(pz-p)(z-z')] . (5.i 6 ) 

J —oo 

If the energy in this fictitious problem is low enough so that the momentum p 2 is below 
threshhold and thus purely imaginary, then it is true that 

F CT (Z) F(Z), (5.17) 

assuming that Z — Z' is "large enough" to kill the exponential. Now, if we take the derivative 
with respect to Z we see that 

F\Z) = f(Z), (5.18) 
F' CT {Z) = 0, (5.19) 

in the low-energy regime where the exponential will be damped because we are below thresh- 
hold. So, we are faced with the very interesting example of two functions which approach 
each other as E — > M + but whose derivatives do not. These CTD effects are relevant for 
the Dirac formalism used here because of the derivative (cr • V) which appears in the lower 
component of the Dirac wavef unctions. This effect is small 

5) The Dirac strengths have a significant real part. The reader may be confused by 
the difference between our approach and the usual procedure of using the optical theorem 
to relate the strength of the optical potential to the cross section. (There the real part 
of the scattering amplitude causes negligible effects.) Here we convert the usual scattering 
amplitudes into their invariant Dirac form to extract the strengths of the scalar and vector 
potentials. Even when the real part of the forward scattering amplitude is neglected, the 
real parts of the Dirac strengths do not vanish. To illustrate this effect, we define a quantity 

Q V (R) = -i [ Z dZ' (F° r - tfj.) p v (B, Z>) 

J —oo v ' 

1 - e ^-p)(z-z')] ( 5 .20) 



X 



as a typical term in the CT equations which will dominate at high energies. Let us further 
define the integrals above to be 
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C R (B,Z)= / dZ' Pv (B,Z') 

J — oo 

x[l-cos(p 2 -p)(Z-Z')], (5.21) 
S R (B,Z) = [ Z dZ'p v (B,Z')sm(p 2 -p)(Z-Z'), (5.22) 

J — oo 

so that 

^(R) = -F^C R (R) - F° r S R (R) 

+iF^S R (R) - iF° r C R (R). (5.23) 

Now, in the previous approach, it was as if F® r were zero. In that case, since C R is positive 
definite, e' Q "' 2 was less than unity. Because our strengths have a significant real part and be- 
cause S R is not positive definite, it is possible to obtain e'^' 2 > 1. Therefore, transparencies 
greater than unity are not disallowed. 

The precise size of this nuclear enhancement effect depends in detail upon the relative 
sizes and phases of C R and S R as well as the strengths of the optical potentials. Thus, 
we cannot predict for certain that the cross section ratios will exceed unity at presently 
feasible energies, although we believe, on the grounds of the above argument, that this must 
eventually happen. 



C. Zero Size 



We also assume that the hard interaction forms a PLC. Then, 

(N\n(q) =Y,r%, m (q)(m\, (5.24) 

m 

= C"*(g)(b = 0| (5.25) 

where r^ m (g) is defined in Eq. QL4|) . Using the specified model, all of the form factors 
are equal and we can neglect the subscripts on Fx, F 2 and on T M (r M = T^^). Then, 
(r M ) _1 r^ m = 1. Later on we shall return to the case where the inelastic form factors are 
not equal to the elastic ones. In particular, in this Dirac case, interesting effects arising from 
a possible different F± and F 2 dependence on the size of the wavepacket are possible. 

It is simple to obtain explicit evaluations of Eq. ( |4.16| ) within the various approximation 
schemes we have described. The results are lengthy; the interested reader can find them 
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in Appendix B. With explicit expressions for the wavefunctions in hand, we use Eq. (|4.15|) 
to calculate the nuclear current matrix element, from which we can construct the nuclear 
tensor and all of the observables. 



VI. DISCUSSION OF NUMERICAL RESULTS 



We have displayed the nuclear current matrix element (NCME), in Eq. ( f4.15|) , as a three 



dimensional integral, with another integral to be performed in the wavefunction itself, for 
each order in perturbation theory. In the basis where d 3 R = BdB dZ d<p, the <ft integral can 



actually be done exactly and analytically [34], but the result is not useful. Therefore, we 



simply perform all of the necessary integrals numerically. We must include 50 integration 
points per integral, using Gaussian quadrature to get stable numerical results. 

The displayed results use kinematics closely related to that of the proposed experi- 
ments [§,0. We do not use precisely those kinematics because our optical potential strengths 
are only known at energies where previous free NN data has been taken. Predictions for 
the upcoming experiments can be readily obtained using simple interpolations. These kine- 
matics are summarized in Table |HI| . For all of our calculations, we take the angle (3 = n/2 
and consider only in-plane scattering. The quantities E\ and 9 e denote the initial electron 
energy and the electron scattering angle, respectively. 

1. Convergence of Exponential Approximation 
From our previous experience in Ref. we learned that the exponential approximation 



(EA) was a good approximation to the full multiple-scattering series. This was seen by 
examining the higher order terms in the various approximation schemes and showing that 
all of the schemes converged to the EA in the energy regime where they were expected to be 
valid. The interpretation of the success of the EA is that the neglect of the path ordering 
includes many more higher-order terms than one would naively expect. We do not expect 
this conclusion to change simply by changing over to a Dirac description. To be careful we 
check to see if this remains true. 
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We calculate the total cross section ratios for the various approximation schemes de- 
scribed above, in the forward direction only to cut down on the computational time required. 
This should be a sufficient test. In Figures § and £| we display the second-order cross section 
calculations in the forward direction, for two values of the internal (quark) oscillator spacing. 
To avoid distraction by the energy dependence of the elementary proton-proton scattering 
amplitude, we fix the optical potential strengths at their Q 2 = 5.96 GeV 2 values. Further, 
we have calculated these results at the values Q 2 = 1, 2, 3, 4, 5, 6, 8, 10, 12, 14, 16, 20, 25 
and 30 GeV 2 . 

We see that the exponential approximation (EA) is a good approximation to the true 
answer. At low energies, the LEE\ follows the EA from the lowest energies to Q 2 ~ 7 GeV 2 
for a light excited state mass and up to Q 2 ~ 12 GeV 2 for a heavy mass. The second-order 
LEE 2 does even better than the LEE\ at higher energies and approaches the EA. It is clear 
that if one were to keep including higher orders in the LEE expansion, we would improve 
the higher energy results. At lower energies, the LEE 2 does not really converge to the EA 
as well as LEE\. While one could argue that the agreement between the EA and the LEE\ 
is accidental, we do not take this view. We believe that at the lower energy end, the LEE 2 
is actually less accurate than the LEE\. This is because we set p/p2m equal to unity in our 
CT potentials, which is accurate at high energies. In a first-order calculation, because of the 
presence of the quadratic terms, we have neglected terms like (p/p^ — 1). In the second-order 
calculation, we have ignored terms like (p/pg — 1)- This kind of effect should be especially 
sizable a low energies. In order to do this calculation more carefully, one should keep these 
"momentum fraction" terms. However, these terms are of the same order as terms which 
have been neglected in the eikonal approximation itself. Thus, to do a really careful job, 
one should also include the first non-eikonal correction to the scattering equations. We do 
not do this here. The only point about the LEE 2 we want to make is that at low energies 
it is approximately equal to the LEE\ and the EA and at higher energies, it more closely 
approximates the EA than does the first-order calculation. 

At the high-energy end, we see that the OBO\ fails terribly, even for l2 C. Of course, we 
expect the OBO to be increasingly worse for heavier nuclei where the propagation lengths 
are longer. However, adding in the second-order contribution to the OBO we find very good 
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agreement for Q 2 > 16 (21) GeV 2 between the OBO2 and the EA, for light (heavy) excited 
state masses. 

At low energies, we really believe the LEE ( the LEE\ in particular) and at high energies 
the OBO is defined to describe the correct physics. What the results in Figures 3 and 4 
tell us is that the approximation schemes overlap so that we can have confidence that the 
EA is the correct answer from Q 2 « 1 to 7 (12) GeV 2 and from Q 2 > 16 (21) GeV 2 for light 
(heavy) excited state masses. Although we have not calculated higher order terms we take 
the position, that with our previous experience in Ref. ||30|| , we can believe that the EA is 
a good approximation to the full multiple-scattering series over the whole energy range. 



2. Integrated Unpolarized Cross Sections 

We next display results for the ratios of integrated unpolarized cross sections. That is, 
we calculate the total cross section, as in Eq. ( |2.19 ), and divide by the Born cross section, 



defined by using a Dirac plane wave for the outgoing proton wavefunction in Eq. ( |4.2j ). See 
Figures 5, 6 and 7. 

The solid curves (circles) are the DWBA cross sections, normalized to the respective 
Born cross sections. The dotdashed curves (diamonds) denote the cross sections, including 
CT effects via the exponential approximation, for oscillator spacings of AM 2 = 1.19 GeV 2 
(M 2 = 1.44 GeV). The dashed curves show the cross section ratios, including CT effects via 
the exponential approximation, for an excited state mass of M 2 = 1.80 GeV. See Table [T 1 1 
for the electron kinematics. 

There are several noteworthy features in the figures. The first is that the optical potential 
strengths, which are calculated from free nucleon-nucleon scattering data, naturally give a 
DWBA cross section ratio which is large at Q 2 ~ 1 GeV 2 and decreases sharply with energy 
before reaching an asymptotic value. This can be understood by noting that the pp cross 
section is only 30 mb at this energy. Thus, we have here an effect similar to the one of 



Frankfurt, Strikman, and Zhalov |62j| . Recall, the NE-18 people took measurements at 
Q 2 =1,3,5 and 6.8 GeV 2 . With an excited state mass parameter of lAAGeV the cross 
section at Q 2 m 3 GeV 2 is the same as at Q 2 « lGeV 2 , although at Q 2 m 5 and 7 GeV 2 the 
cross section rises rapidly with energy. 
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We want to stress, at this point, that the precise value of the oscillator parameter, which 
controls the value of the excited state masses in the CT models, is unknown except for the 
fact that is on the order of hadronic mass differences. 

Therefore, by simply increasing the value of the first excited state mass, we can postpone 
the onset of transparency to higher energies, see the dashed curves. The suppression of the 
transparency for an excited state mass of 1.80 GeV 2 may be too great to agree with the SLAC 
data at Q 2 « 3 GeV 2 . Thus, it seems that an oscillator spacing such that M 2 ~ 1.6 GeV 
would give a cross section ratio which would appear to be Q 2 -independent if one only looked 
at Q 2 = 1,3, 5, 6.8 GeV 2 . 

The SLAC data are taken at kinematics slightly different than what we use. However, 
for cross section ratios, the differences in the electron kinematics lead to essentially the same 
results. 

That the SLAC experiment has seen CT is certainly possible, given the energy depen- 
dence of the elementary pp cross section and the optical potential strengths. However, before 
one can be absolutely sure, it is desirable to take more data points and see what happens 
at, say, Q 2 ~ 2 GeV 2 . It would also be nice to increase Q 2 and really see a dramatic rise in 
the ratio a/a B . 



3. Integrated Longitudinal and Transverse Responses 

It has long been known that the longitudinal response function is suppressed relative to 
the transverse response in the quasi-elastic region pB|-p^]. The theoretical interpretation of 
this quenching has been difficult and no consensus, beyond the failure of naive applications 
of traditional nuclear theory, has yet been achieved. In particular, it has been suggested 
that a modification of the properties of the nucleon in medium or the presence of non- 
nucleon degrees of freedom is responsible for this suppression. It has also been suggested 
that relativistic effects may play a role. Indeed, Do Dang and Van Gai |kJ have shown that 
the presence of large scalar and vector potentials in the nucleus may be responsible for this 
effect. Since our calculation is based on such a scalar- vector phenomenology, it is legitimate 
to ask whether or not our calculation displays this quenching. 

As the integrated cross sections only involve the integrated longitudinal and transverse 
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response functions, Eq. fl2.19|) , it is useful to examine the ratio of 



K L = Jd( sm(R L (Q _ [dq± q± g^(g±) ,„ ^ 

TZ T J d( sin C Rt(C) ~ / dq ± q± R T (q±) ' 1 ' ] 

where ( is the angle between the virtual photon and the ejected proton and q± = |p — q| ~ 
|q|sin£. Since the response functions Rl and Rt go to zero rapidly for q± > 1.5 /to -1 , 
we take the limits of integration in q± to run from zero to some number q™ ax pa 2 /to -1 . 
The contribution to the angle integrated responses TZl(t) from gj_ greater than 2 /m -1 is 
negligible. Note that this integration is over angle and not over the energy transfer. Thus, 
our ratio TZl/IZt has very little to do with the integral involved in the Coulomb sum rule. 

If longitudinal quenching occurs, IZl/TZ-t < 1- The dependence of this ratio as a function 
of Q 2 , for quasi-elastic kinematics, can be inferred from the concept of y— scaling [|i6| ,|6"9|-|75 



In order to determine the Q 2 dependence, we define the new reduced response functions 
and tt as 

n L (q, g ) = ^Gl(Q 2 ) (jf\ rz(q, 9b), (6-2) 



A ^2 ( Q 2 _ 
r2 



^T(q,go) = ^Gl^Q 2 ) ) r T (q,g ), (0.3) 



where Ge and Gm are the electric and magnetic form factors of the nucleus, which are 
weighted sums of the electric and magnetic form factors of the constituent protons and 
neutrons. The statement of y-scaling is that at large enough momentum transfers, the 
reduced response functions ri, and tt are functions only of the variable 



|q| . M N q Q 

y = — tt + 



The quantity y is to be interpreted as the longitudinal (parallel to q) component of the initial 
nucleon momentum (in Born approximation) when the transverse (to q) component is zero. 
We note that many choices of the scaling variable exist. At the top of the quasi-elastic peak 
(x = Q 2 /(2Mjvg°) = 1), y = for all values of Q 2 . Thus, in this kinematical limit, which 
is precisely where we have calculated, the functions and are independent of Q 2 (when 
y-scaling is valid). In principle, ri(y) and rr{y) should be equal. Experiments have found 



that although tl and rr seem to functions of y only, they are not the same function [73 
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Thus, simply from looking at the remaining Q 2 dependence in Eqs. ( |6.2|) and ( |6.3|) , we see 
that 

— - constant, (6.5) 

from above. This is indeed what we see in our calculation. 

In Figures 8, 9 and 10 we display the ratio IZl/TIt as a function of Q 2 for 12 C, 40 Ca 
and 208 Pb. As before the circles, joined by solid lines, show this ratio for the DWBA 
calculation. The dotdashed curves (diamonds) and dashed (boxes) curves show this ratio 
for the CT-included cases (via the EA) for the light mass (M2 = 1.44 GeV) and the heavy 
mass (M 2 = 1.80 GeV) cases. We also display the result for the plane wave or Born case, 
which is shown in the figures as the fancy boxes joined by the dotted curves. 

At higher energies the CT-included ratios are slightly more suppressed than the DWBA. 
But only slightly. Thus, even if a careful Rosenbluth-type separation could be performed, 
the differences in TZl/TZt predicted by the Glauber and CT calculations are too small to 
experimentally differentiate between the two. 

4- Gauge Invariance and Current Conservation 

We have mentioned earlier that our calculation violates gauge invariance and current 
conservation (CC). One of the main causes is that, in our NCME, we treat the bound state 
and scattered state as derived from different Hamiltonia and thus are not members of the 
same complete set. This is a typical problem for distorted wave calculations of this type. 
(For the initial- and final-state wavefunctions to be derived from the same Hamiltonian is 
necessary but not sufficient to obtain gauge invariance.) Gauge invariance is often imposed 
artificially fl33l , |76l . In our view this is not the best solution, since one must be very careful 



to not introduce unphysical terms into the calculation. It is true that for our calculation to 
make any sense at all, the current must be at least approximately conserved. This subsection 
examines the question of CC. 

There are many ways that one can consider quantifying the amount of current conserva- 
tion violation (CCV). The nuclear current has four complex components, so directly looking 
at q^J^ seems tedious. Instead, we consider the relevant elements of the nuclear tensor, 
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W 1 *" . We would like to investigate how the CCV's change as a function of Q 2 . Since W^ u 
is a function of angle (q±), we define 

W V = ^ [dq ±q± W" u . (6.6) 



iqi 

where the energy- integrated nuclear tensor W^ u is defined by Eq. ( |2.5|) . The quantity q± = 
|p— q| « |q| sin £, where £ is the angle between the the virtual photon and the ejected proton. 
Choosing the Z-direction to be defined in the direction of q, gauge invariance requires 

AR L = %V _ 00 g3 _ 33 (6.7) 

to be zero. We display results for AR^ as a function of Q 2 for our three canonical nuclei, 
12 C, 40 Ca and 208 Pb, and for both internal oscillator spacings, M 2 = lAAGeV (light) and 
M 2 = 1.80 GeV (heavy). 



We see, from looking at Figures fTJ, [12] and |13| that our calculation violates current 
conservation at the 10% level for the DWBA at the lowest energies but only at the 3-4% 
level for the CT cases. At higher energies, our calculation becomes more and more current 
conserving. Notice that the current is mostly conserved for the plane wave cases, becoming 
less so for the heavier nuclei. 

Of course, in nature, current is conserved exactly. However, we make the provocative 
statement that the improvement of current conservation in our calculation over the tradi- 
tional Glauber treatment is more evidence that CT effects are necessary for an accurate 
theoretical treatment of (e, e'p) reactions. 



5. Differential Cross Sections and Normal Polarizations 

In this section, we display the results from our calculations for the differential cross 
sections (as a function of the transverse momentum) and of the normal polarizations. These 
are the quantities which will be measured in the upcoming CEBAF experiment We 
note that these are the first predictions of these observables at these energies which include 
the effects of CT. It is important to recall that, in the CT limit, the normal polarization 
should vanish. 
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We only display here, for reasons of space and continuity, some representative results for 
12 C nucleus. The results for the other energies which we calculate and for 40 Ca and 208 Pb 



are presented in Ref. ]77|]. Here, the quantity Ei is the initial energy of the incident electron. 
Different values of Ei affect only the kinematical weights used to construct the cross section 
and polarization; the nuclear response functions are independent of the electron kinematics. 

Clearly visible, in Figures 14 and 15, in the shapes of the differential unpolarized cross 
sections is the shell structure of the nucleus. For instance, 12 C has 4 p-shell protons and 
only 2 s-shell ones. The wavefunctions for the s-shell nucleons peak at the origin, of course, 
while the p-shell ones have a node there. This is the reason why the cross section for 12 C 
has a maximum at about q± ~ 0.5 ]mT x . Similarly, the shell structure for 40 Ca and 208 Pb 
can be discerned from the figures in Ref. |77j . 

We predict that the energies proposed in the experiments are not high enough to see 
the normal polarization vanish. However, in 12 C there is a measurable suppression of the 
polarization for the CT-included case of light excited state mass (M 2 = lAAGeV). A 
detailed comparison with the results of the SLAC NE-18 measurement of cross section 
ratios is necessary to see if using the light excited state mass is still viable. Our preliminary 
results are that using such a low mass may not be inconsistent with the SLAC NE-18 data. 
In m Ca and 208 Pb, there is only a small suppression of the normal polarization even at 
Q 2 = 20 GeV 2 . 

Thus, it seems that the normal polarization is a difficult quantity in which to observe 
CT effects. But a precise experiment could be successful. 

Note that our formalism does indeed yield vanishing normal polarization at high energies. 
To see this we look at the observables at Q 2 ~ 300 GeV 2 , with optical potential strengths 
fixed at their Q 2 = 5.96 GeV 2 values. At this extremely high energy, the EA wavefunction 
is essentially the same as the OBO\ wavefunction and the normal polarization is seen to be 
essentially zero for l2 C over the range in q± we consider, although for m Ca and 208 Pb the 
polarization is still sizeable at the larger angles, even at this enormous energy. 

In Figure 16 we show the predictions from the DWBA and the Born approximation 
(BA). The Born cross section is larger in magnitude than the DWBA while the polarization 
is identically zero for the BA. The dotdashed curve is the prediction from the exponential 
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approximation, as before. The dotted curve is the prediction for the LEE\ and the dashed 
curve is for the OBO\. In all the CT curves, we have used an excited state mass parameter of 
M 2 = 1. 44 GeV. Much of the differences between the figures for different electron kinematics 
is due to the huge difference in the Mott cross sections, see Table |T|. There is also a sizeable 
difference in Vr- However, it is important to note that all of the Q 2 = 5.96 GeV 2 results lead 
to cross section ratios which are the same to better than 1%. Since absolute magnitudes of 
the cross sections are larger for the greater incident electron energies, due to the Mott cross 
section, these kinematics are preferred experimentally. It is also interesting to comment 
on the OBO\ curve. We see that the OBO\ is, already at Q 2 = 5.96 GeV 2 , predicting a 
suppression of the normal polarization. We know that the OBO\ is definitely inaccurate at 
these low energies. However, the high-energy nature of the approximation manifests itself 
in giving small polarizations. 

6. Differential Unpolarized and Normal Response Functions 

In this section we display representative graphs of the separated response functions, 
see Figures 17 and 18. Using these response functions, the total differential cross sections 
and normal polarizations can be constructed by summing them together with appropriate 
kinematical weights, given by the Vs. 

Since these quantities are more difficult to extract than the total differential cross sections 
and polarization, and Saha MM is not planning to do so, we only display a few representative 
figures. The reader who is interested in more results of this kind are referred to Ref. |77| . 

We see from looking at the response functions that nothing really special happens for 
the "unpolarized response functions" . Response functions represent spin-dependent or spin- 
independent measurements. It is convenient to adopt the phrases "unpolarized response" 
and "polarized response" or "normal response" to indicate the spin-independent or spin- 
dependent (dependent on normal component) responses. 

In lead, there is strong absorption and the effects of CT do not manifest themselves 
until very high energies, making the experimental verification of CT in lead unlikely if one 
looks only at unpolarized observables. However, in the normal responses, we see something 
quite different and quite interesting. In the transverse normal response, Klf, we see that for 
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both light and heavy excited state mass cases there is a huge enhancement in the region of 
q± « 0.4 — 1.5 fmT 1 . That is, a ratio of CT response to DWBA response ranges all the way 
from zero to infinity over this angular range. 

We note that this CT effect in Kj, is not of the asymptotic nature one normally thinks 
of as CT, lack of absorption, something going to zero, etc. However, if one takes the CT 
wavefunctions we have derived seriously, then the low energy spin observables are affected 
in this way. In this sense, this enhancement in K!f is a pre-asymptotic effect of CT. This 
effect is also present in 40 Ca and in l2 C but to a lesser, and still lesser, extent. 

If detected, this enhancement in BJf would be an unambiguous signature of CT. We 
stress that this enhancement begins at momentum transfers as low as Q 2 ~ 2GeV 2 (!) and 
continues to the highest energies. However, the separation of this response function from the 
total normal polarization may be difficult. The kinematical weights, at these energies and 
angles, are all approximately equal, see Table [IT], so that the transverse normal response 
only contributes about one part in twenty to the polarization. The predictions presented 
here indicate that attempting this separation may be worthwhile. 



In all of the above calculations we have assumed that the wavepacket formed in the initial 
hard collision was of exactly zero size. This is not realistic. We do expect, however, that the 
size of the wavepacket will decrease with increasing momentum transfer. It is interesting 
to take this non-zero size into account. This is done by first considering the form factors 



part factored out. Further, we assume that the Dirac and Pauli form factors have the same 
internal dependence. That is, we take 



A. Other Applications 



1. Finite Size Effects 



defined in Eqs. ( f4.3|) and ( |4.4| ) to be the product of the usual form factors with the internal 




(6.8) 



(6.9) 
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In the absence of a compelling reason, note that we choose the internal dependence for F% 
and F 2 to be the same. With this simplification, the internal form factors simply factor out 
and (no sum over //) 

[^(g)]- 1 r^ 2m (g) = ^^, (6.10) 

jN,N{q ) 

which is proportional to the unit matrix in Dirac space. Now, to first order in the interac- 
tions, we identify as the relevant quantities 



7/(2i)/ /-R 7 7 I\ _ \ - fN,2m(q 2 ) 

U cso(so)\ a i Z > Z ) ~ 1^ — f 7~-K 

m JN,N{q ) 

m,2j 



xfflB.^M (6.11) 



where it is understood that, to the present order, only m = 0,1,2 are nonzero. In order 
to obtain explicit expressions for the internal form factors, we need to further specify the 
internal part of the hard scattering operator. It is generally believed that the size of the 
wavepacket is inversely proportional to the momentum transfer of the virtual photon. Sup- 
pose a nucleon, comprised of three quarks absorbs a photon of three-momentum q. Each 
quark, then, recieves approximately |q|/3 of the photon's momentum. So, we use a form |78 
for the internal part of the hard operator 

7%{q) = T£(g)e^ l/3 . (6.12) 

This form is suggested by caricatures of perturbative QCD calculations; see Eq. (1) of 



Ref. [24 



In that case, the internal form factors /n,n, /jv,2 and /jv,4 are given by 

fN,N( x ) — 1 ~ xV^e* 2 erfc(x), (6.13) 
f N>2 (x) = -x 2 + (x 3 + ~x) y/ne x2 eTk{x), (6.14) 

e 1 \ 1 4 1 2 

JnA x ) = 2 X + 1 X 

+ (-\x b - ^x 3 + -x) ^fKe x2 eiic(x), (6.15) 
V 2 2 8 / 

where x = x(q 2 ) = |q|6#/6. 

The results for non-zero initial wavepacket size are displayed in Figures 19 and 20. In 
Figure 5 we show the ratio of cross sections for 12 C as a function of Q 2 . As before, the 
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solid line denotes the DWBA (where finite size of wavepacket has no effect) and CT is 
included via the exponential approximation for both light (M 2 = 1.44 GeV) and heavy 
(M 2 = 1.80 GeV) excited state masses. We have also overlayed the zero size calculations for 
comparison. For both the light (diamonds, dotdashed) at heavy mass (boxes, dashes) cases, 
the curves including the finite size are slightly below the zero size cases. That the effect is 
small even for the case of the light excited state mass is an indication that assuming exactly 
zero size is not a bad approximation. 

In Figure 20 we show the CT curves for the current conservation violation (CCV) quantity 
ARl as a function of Q 2 . We label the curves E A to indicate that the CT effects are 
calculated in the exponential approximation. The superscripts LT and HVY indicate the 
mass of the excited state, light (LT) or heavy (HVY). The subscripts FS stand for "finite 
size" and those curves have included these finite size effects in the manner described above. 
Thus, the dotdashed (dashed) curves should be compared with one another. Looking at 
Figure 20 we see that the effects of including finite size improve the situation as far as 
gauge invariance is concerned. Including the size of the initial wavepacket shows that the 
CT calculations are gauge invariant at the 2% level at low energies whereas the zero size 
calculations are current conserving at the 3% level, at low energies. However, the trend to 
make the calculation better obey the requirements of current conservation is encouraging. 
We can conclude that assuming the initial wavepacket to be of exactly zero size introduces 
only small errors and can therefore be considered a good approximation. 

We do not display results for the finite-sized calculation of TZl/TZt because there is no 
discernable difference in this quantity between the finite-sized curves and the zero-sized ones. 



2. Fermi Motion 



It has recently been argued [|^[79H 81|| that the Fermi motion of the nucleons bound 
in the nucleus strongly enhances the effects of color transparency. Suppose, in the (e, e'p) 
reaction, that the virtual photon three-momentum is labelled by q and the detected proton 
three-momentum is called p. Suppose that we call p = q + k. Sometimes k is called the 
"momentum of the struck nucleon". This is only true in Born approximation (the CT limit) 
because final-state interactions can influence the outgoing proton's longitudinal momentum. 
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In exactly quasi-elastic kinematics, the bound nucleon is treated as at rest and free and so 
k = 0. Fermi motion allows (requires) k ^ 0. In the above calculations, we have assumed 
that k = qj_ where is very small in magnitude, |p| ~ |q|, and is purely transverse to q. 
This assumption leads to the differential distributions in Section [Vl| and in Ref . |77] . 

It has been shown [fTS|-|ST|] that the component of k which is parallel (or anti-parallel) to 
q has a huge numerical effect on the calculated transparencies. This section is devoted to 
an investigation of this effect in our models. 

We let the component of k which is in the q direction to be called fen Thus, k = fenZ-|-q_i_. 
When q^ = then fen = y, where y is the scaling variable introduced in Section [VI 3| . Non- 
zero values of fen are especially important in (p,pp) reactions [^2j because the experimental 
results of Refs. []3],^| are presented in terms of bins of fey. As noted earlier, a correct ac- 
counting of this effect, in the kinematics of the Brookhaven experiment, leads to a better 



understanding of the data [80]. 

We have calculated the cross sections and other observables for Q 2 = 0.96 and 20.86 GeV 2 
for values of fen = —150,-75,0,75 and 150 MeV. In particular, we assume that p • Z = 
|q| +k\\, where q is in the Z— direction. The results for the total cross section ratios, TZl/TZt, 
and the CCV measure ARl are shown in Figure 21. In all of the figures we display four 
curves. The dotted curve is the Born approximation, the solid curve is the DWBA, and the 
dotdashed and the dashed lines are CT results for light and heavy excited state masses. In 
the plots for the cross section ratios, we display the ratios for the cross section at fey divided 
by the Born cross section at that same value of fey, except for the Born curve, which is the 
cross section at fey divided by the Born cross section at k\\ = 0. At low energies, the effect of 
finite fey is pretty small on the cross section ratios, about 20%. At high energies, the effect is 
quite large for the CT cases, although still a small effect for the DWBA. However, nonzero 
values for fey lead to a large reduction, over a factor of 2, in the Born cross section at low 
and high energies. 

Another extremely interesting result from these calculations is shown in the measure 
of CCV. In particular, we see that as fey moves away from zero, current conservation is 
violated more and more, reaching the 40% level for |feJ = 150 Me V, at low energies. At 
higher energies, the situation improves so that the violations are at the 5% level. Clearly, 
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violations of 40% or more are intolerable and our calculation probably cannot be trusted 
there; we shall come back to this point later. Violations of 5% probably do not affect the 
physics too much. Thus, there is no problem in trusting our calculation at high energies or 
at small values of k\\. 

What is the cause of the big change in CCV, due to the non-zero fcy? One can get a hint 
by looking at the nuclear current for scalar electrodynamics. In that case, the current can 
be written as 

J^oc j d'p^ip + pTMp'), (6-16) 

where <f>i(p) is the momentum space bound state wavefunction and 0/(jo') is the scattered 
wavefunction. In the Born approximation where <f>f(p') = S 4 (p' — p — q), the integral can be 
done simply. The initial nucleon momentum p can now be indentified with k, defined above. 
Then, the statement of current conservation is that 

q^cxk\\\q\, (6.17) 

should vanish. We have used that k^ w M% and chosen x = Q 2 / (2M N q°) = 1. There- 
fore, we see that for vanishing k\\ the current should be exactly conserved. Including the 
spin degrees of freedom makes this no longer exactly true, but this explains why the Born 
approximation is the most current conserving of our models. This also explains why, as 
fcy gets bigger, the violations of current conservation also grow. For large Q 2 , we find that 
g M J M oc Q 2 k\\/M N [1 + 0(Mff/Q2)}, which explains why the calculations become more gauge 
invariant at larger energies. For non-zero ku the above argument indicates that the current 
is not conserved. Why not? Obviously, something is missing in our calculation. The most 
important cause of CCV in our calculation is probably coming from the lack of orthogonality 
between the initial state and final-state wavef unctions. Indeed, we have left out effects in 
the initial state, such as particle-hole excitations as well as explicit isobar degrees of freedom 
which can alter the nucleon-nucleon force. Also, we have ignored the possibility that the 
nucleons are altered at all in the medium. Further, our assumption that the current operator 
is the free nucleon operator is probably not warranted; off-shell effects may be important. 
A more careful treatment of some of these effects could lead to improvements in the gauge 
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invariance of our models. The source of these CCV's are of much interest and deserve further 
study in the future. 



3. Lower Components 

We examine the effects of the lower components of our model wavefunctions here. We 
only present results for the light mass case and for the 12 C nucleus. 

First, we artificially turn off the lower components in the scattered waves. The results for 
the integrated cross section ratios are shown in Figure 22. The solid and dotdashed curves 
are calculated using the the DWBA and EA wavefunctions with the lower components 
turned off, divided by the Born calculation also with no lower components. The dotted 
curve in Figure 22 is the ratio of the Born calculation with no lower components to the Born 
calculation with lower components. The figure shows a striking effect: the lower components 
are not significant in the predicted cross section ratios. These cross section ratios are seen to 
be a quantity which is very insensitive to the details of the wavefunctions. This is actually 
quite amazing since the Born calculation changes so much when the lower components are 
turned off. 

However, we point out that these lower components are crucial to the approximate 



satisfaction of current conservation, see Figure |23|. Further, turning off the lower components 



has a huge effect on the integrated ratio of longitudinal to transverse response, see Figure |24. 
Thus, we see that the suppression of the longitudinal response relative to the transverse 
response is really entirely located in the lower components. This is consistent with the 
notion, mentioned earlier, that this longitudinal quenching was due to relativistic effects. 
Naturally, since this ratio is quite different than previously, all of the other observables such 
as differential unpolarized cross sections, differential normal polarizations and differential 
response functions are also altered by a large amount. 

We have also performed a calculation in which we take the ratio of lower to upper 
components in the scattered wave to be the same as for a free Dirac plane wave, %? M . This 
replacement has essentially no effect on any of the observables at the energies we consider. 
Only the ARl is altered noticeably, with the DWBA violating the current conservation a 
little more (~ 5 — 10%) without lower components than with them and the EA conserving 
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the current a little bit more. These results show that the CT derivative effects, described in 
Section |V B , are not numerically significant and practically all of the physics is independent 
of these CTD's in the lower components. Since the differences are small we do not display 
any numerical results for this calculation. 

We also investigate turning off the lower components in the bound state wavefunction. 
This change introduces only very small corrections to the observables and has only a small (< 
5%) effect on the level of current conservation violations; the changes in the other observables 
are not important. This is consistent with our expectation; the QHD bound state lower 
components contribute only 2 — 3% to the total probability. That is, / dr\F\ 2 / J dr\G\ 2 ~ 
0.02 — 0.03. We do not display any numerical results for this calculation, either. 

We conclude that the lower components are essential to our description of the ejectile 
wavefunction. Without these inherently relativistic components we would lose approximate 
current conservation and the experimentally observed quenching of the longitudinal response, 
although the integrated cross section ratios remain unchanged. However, examination of the 
details shows that the dominant effect in the lower components of the wavefunctions is the 
a ■ p term and that the complicated terms which appear in the explicit evaluations of the 
wavefunctions are unimportant. 

We have also explored the effects of changing our bound state wavefunction. Instead 



of using the QHD wavefunction of Horowitz and Serot ||38|| , we use harmonic oscillator 
wavefunctions. Using these bound state wavefunctions leads to only very small changes in 
the cross section ratios and in the violations of current conservation. However, the differential 
cross sections are slightly different than those presented in Section [VI 5| . This is because 
the p— wave wavefunction has a smaller (R 2 ) (higher momentum components) than the 
corresponding QHD wavefunction. As a result, the differential cross section calculated from 
this oscillator wavefunction peaks at a larger value of q± and has smaller contribution at low 
q±. The contribution to the differential cross section from the s— wave fall faster than the 
increasing contribution from the p— wave, so that the bump that we saw in the QHD case 
disappears. 

It is possible to fix the (R 2 ) for the harmonic oscillator wavefunctions to agree with that 
predicted by the QHD model. Doing this improves the agreement between the p— waves but 
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at the expense of the agreement between the s— waves. Thus, the differential cross sections 
seem to be sensitive to the precise form of the bound state wavefunctions. 

However, the main conclusion of this work on polarization, that the normal polarization 
is not hugely suppressed at present energies, remains unchanged when harmonic oscillator 
bound states are used. 

VII. SUMMARY AND CONCLUSIONS 

The effects of color transparency (CT) in quasi-elastic (e, e'p) reactions at large mo- 
mentum transfers are explored here. We have focused on the (e, e'p) reaction since lepton 
scattering is simpler than hadronic scattering. Of course, CT effects should also be present 
in other reactions such as (p,pp), (ir, ftp), etc. The methods presented in this paper can be 
applied to these reactions. 

In this paper we include the effects of proton and photon spin. Thus, we call this work 
"Vector CT", reflecting the vector nature of the photon. So far no published work has 
included the effects of spin in calculations of CT. The main motivation for undertaking 
this calculation is the proposal by Saha and collaborators |||| who plan to measure the 
normal component of the proton polarization in (e, e'p) experiments. This polarization is an 
interesting quantity because it vanishes in the absence of final-state interactions. Further, 
the photon really is a vector particle and its spin should play some role in the scattering. 
Dirac phenomenology is used to construct the DWBA. The internal operators and states are 
then embedded in to this formalism. We obtain detailed expressions for the CT wavefunction 
in terms of these internal operators. 

In order to calculate the CT effects, it is necessary in this approach to choose an explicit 
form for the wavepacket-nucleon interaction. Although this interaction is not precisely 
known, some general properties are generally agreed upon. We choose a simple representative 
form for this interaction which is consistent with the known constraints. 

Further, it is also necessary to assume a model for the internal space on which the CT 
operators act. For simplicity, we choose the baryon spectrum to be represented by a two- 
dimensional transverse harmonic oscillator. With this choice comes a single free parameter, 
the oscillator spacing, which determines the masses of the nucleon resonances. This spacing 
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is characterized by the mass of the first even parity excited state, which we label by Mjv*. 
In this paper, we have chosen the values M N * = lAAGeV and M N * = 1.80 GeV for the 
numerical results we display. Increasing this mass postpones the onset of CT to higher 
energies. 

We evaluate the formal distorted wavefunction by making various approximations to 
sum up different parts of the multiple-scattering series. In each energy regime we define 
an unperturbed piece of the path-ordered exponential, which we can solve exactly, and a 
perturbed piece, which we treat systematically in perturbation theory. We have developed 
high-energy and low-energy approximations. However, the different approximation schemes 
overlap so that we are confident that we can accurately approximate the exact distorted 
wavefunction at all energies (above 1 GeV). In particular, we find that simply ignoring 
the path ordering symbol in the wavefunction serves as a good approximation of the full 
distorted wavefunction for all energies. 

Although we use over-simplified models for the wavepacket-nucleon interaction and the 
internal baryon space, the methods of calculation presented in this paper are more general 
and can be used with more realistic interactions and models. 

We calculate ratio of cross sections for 12 C, 40 C a and 208 Pb as a function of Q 2 . In 
contrast to earlier calculations which assume a constant value of 40 mb for the elementary 
proton-proton cross section, we take the cross section and other optical potential strengths 
directly from data. Thus, we confirm the assertion of Frankfurt, Strikman and Zhalov, that 
energy dependence in a/a B is expected. This is because the pp cross section, in the energies 
of the experiment varies in such a way that the Glauber treatment decreases the transparency 
ratio. Therefore, since the preliminary results of the experiment, so far, see a small variation 
in the ratio as a function of Q 2 , there may be some transparency effects to compensate. We 
conclude that the excited state mass which is consistent with the preliminary data is about 
Mjy* 1.6 GeV or more. 

A typical problem with calculations of this sort is the lack of current conservation. We 
show that the violations of gauge invariance are at the 10% level for the DWBA and only 
at the 4% level for the CT wavefunctions at low energies but are all at the 1 — 2% level 
at the highest energies, Q 2 ss 20 GeV 2 . Inclusion of CT effects is therefore desirable, even 
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at low energies, because it improves the situation from the DWBA with respect to current 
conservation and gauge invariance. 

Because of the extra spin degrees of freedom, there are many many other quantities which 
we calculate. We display results for the ratio of integrated longitudinal to transverse response 
as well as for differential cross sections and normal polarizations. The differential cross 
sections do display a noticeable, and detectable, increase at large momentum transfers. This 
is just the data which is summarized in the integrated cross section ratios described above. 
In the normal polarizations, we see that the deviations from the usual Glauber treatment do 
not seem to be significant or measurable at any energy. Although Q 2 20 GeV 2 we see, in 
12 C a moderate decrease of the polarization for the light mass CT case but not for the heavy 
case. We conclude that the energies at which the normal polarization should completely 
vanish are, unfortunately, quite high. Frankfurt, Strikman and Zhalov have obtained a 
similar result (private communication). 

We have also calculated the individual separated response functions, eight in all for in- 
plane scattering, which require many measurements at many angles in order to accurately 
extract from the cross sections and polarizations. We notice, in some of the response func- 
tions, that some pre-asymptotic behavior of the CT wavefunctions are manifest are possible 
detectable. This is an effect which increases in heavy nuclei like 208 Pb. In particular, if it 
is possible to separate out from the normal polarization the responses Rjj-, Rti Rtt ^^-d 
R^ T , then it should be possible to see in R% a huge enhancement over the Glauber result at 
|q_i_| ~ 0.5 — 1 jmr x for Q 2 > 2GeV 2 \. This extraction may be difficult because, although 
the kinematical weights are all approximately equal, R% contributes only about one part in 
ten at most to the total polarization. 

Further, we have also investigated the effects of non-zero wavepacket size and of Fermi 
motion. Including the non-zero size of the initial wavepacket produces only small changes so 
its neglect is a very good approximation. This is because of effects of wave packet expansion. 
Very small sized wave packets expand quickly into small sized wave packets. However, a 
slight improvement in the current conservation is obtained by including the non-zero size. 
The Fermi motion, on the other hand, can have a big effect on the cross section ratios if the 
"initial nucleon momentum" is along the direction of q. However, for values of k\\ > 150 
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MeV/c, which is this component, the violations of current conservation become severe at 
low energies and our calculation cannot be trusted there. At high energies, the calculation 
becomes more gauge invariant and we see an interesting effect for non-zero k\\. This effect 



is especially important for the (p,pp) reaction where the data is displayed in bins of k\\ |80 

Lastly, we have also investigated the effects of the lower components on cross section 
ratios, current conservation, polarization, etc. The complicated nature of the DWBA and 
CT lower components are not significant and taking the ratio of lower to upper components 
for these distorted waves to be the same as for plane waves introduces no discernable results 
(CTD effects are unimportant). The bound state lower components are similarly seen to be 
unimportant in the prediction of most of the observables. Finally, we have investigated the 
effects of using harmonic oscillator bound states instead of the ones given by QHD. The cross 
section ratios are essentially the same as for the QHD bound states, while the differential 
cross sections, polarizations and response functions all differ somewhat from those using 
QHD bound states. 

In conclusion, we can make the following short comments. The recent SLAC data im- 
poses constraints on the allowable models shown in this paper. Despite the lack of significant 
Q 2 variation in the data it is possible, because of the energy dependence of the elementary 
pp observables, that the data actually may be an example of the manifestation of CT. A 
measurement of the normal polarization in (e, e'p) reactions does not seem to be a good 
way to see CT effects at moderate Q 2 . However, a measurement of the normal transverse 
response in a heavy nucleus such as 208 Pb does seem to afford the opportunity to see CT, 
unambiguously, at quite low momentum transfers. 
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APPENDIX A: EXTRACTION OF OPTICAL POTENTIAL STRENGTHS 

We describe how we obtain the strengths of the scalar and vector optical potentials, first 
shown in Eq. (|3.9|) , and summarized numerically in Table |I[ In particular, we extract these 
strengths from proton-proton elastic scattering data. In order to show exactly what we have 
done, it is useful to recall some aspects of nucleon-nucleon scattering. 



1. Nucleon-Nucleon Scattering 



We use the parameterization of the nucleon-nucleon scattering amplitude of [f47H55 



= A + B(Ti ■ <J 2 + iqC {(Tin + cr 2 n) 



+D<Ti ■ qcr 2 • q + Eai Z cr 2Z . (Al) 



The amplitudes A, B, C, D, E can be easily related to the s-channel helicity amplitudes |82 



One can also parameterize the NN scattering amplitude in a Lorentz invariant 
form [f47H53|], by describing the spin dependence in terms of Dirac matrices, 



F = F S + F vl Zl2» + F t a? u a 2flu + F p7l 5 72 5 

+F al hh^ (A2) 

Naturally, the quantities F s , F v , etc. are functions of the momentum transfer, q. By 
taking the matrix element of this operator and doing the Dirac algebra, leaving F to be 
an operator in the Pauli spin space, we can equate the two descriptions, modulo some 
kinematical factors. These factors are known and are derivable from the relativistic impulse 
approximation fTFH53f. In fact, the relationship is, in the center-of-mass frame, 

t t f c 
X K X K 2ik XAlXA2 ' ~ 

u{k', X[)u{-k', A' 2 ) F u{k, Ai)«(-k, A 2 ) (A3) 
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where \\ are Pauli spinors for helicity A and w(k, A) are free Dirac spinors corresponding to 
momentum k and helicity A; k is the center-of-mass momentum. 



2. The Lorentz Invariant Forward Amplitudes 

We are concerned with the case of elastic scattering of protons from nuclei. It is in 
this case that we can extract the optical potential strengths. Within the context of the 
relativistic impulse approximation, the T-matrix for nucleon-nucleus scattering is simply a 
sum of the t-matrices of the elementary nucleon-nucleon scattering. If we neglect nuclear- 
medium modification to the bound nucleons, off-shell effects, and 1/A corrections, the optical 



potential is approximately equal to the impulse approximation T-matrix [ 49(1 . In a spin 
saturated nucleus, the matrix element of this Dirac operator involves a trace over struck 
nucleon spins which eliminates all terms except the tensor, scalar and the time component 
of the vector terms. Thus, in this relativistic impulse approximation, the optical potential 
is of the form, 

a ' q -F t (q) Pt (q)}, (A4) 



M N 

where p s , p v and p t are scalar, vector and tensor nuclear densities and pi a b is the lab frame 
momentum. The kinematical factors on the right hand side of the above equation, r = 
—Airipiab/MN, come from the relationship between the Lorentz invariant Feynman amplitude 
and the usual description of the optical potential in the impulse approximation. 

The coordinate space optical potential is obtained by Fourier transforming the momen- 
tum space potential. In the limit of very high energies, the scattering is predominantly in 
the forward direction. Thus, we can approximate F s (q) and F v (q) by their forward direction 
values Fg and F®, and neglect the tensor term since q is small. Thus, we arrive at f47| , |4"9| 



U opt (R) 



^p 8 (R)+^ Pv (R) . (A5) 



In the forward scattering approximation, the relationship between the Pauli and Dirac 
amplitudes can be written as a matrix equation | 48fl . Since A, B,C, D, E are to be taken 
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directly from data, we can simply invert this matrix to obtain the scalar and vector density 
strengths from the Pauli amplitudes. In particular, we find 

' f £(2£ + 1)A° + (2£ 2 - 1)B° 

(A6) 



2^(1 + 



+2M 7V £(2£ 2 - 1),/P|C + ^r-^E 



Z - 1 2(f 



" 2e 2 (i+o 



-B°- 2M N d^C° 



(A7) 

where £ = ^fsjlM^ is dimensionless Notice that the quantity D° does not enter in for the 
scalar and vector strengths. It is interesting to realize that inasmuch as the central potential 
dominates, it is really only the quantities A and C° which are important. Indeed, if the 
differences between the scalar and vector nuclear densities are neglected, a relevant quantity 
which appears in the Dirac eikonal formalism is 

2^p^a6 

That this should be true is not a surprise. A is just the forward scattering amplitude, whose 
imaginary part is proportional to the total cross section. 

3. Extracting Forward Amplitudes from Data 

The problem is how to extract the NN parameters A°,B°,C°, E° from the data. We do 
this in two ways. In the first way, we take some recent phase shift solutions to reconstruct 
the helicity amplitudes and then to use the above relations to get at the quantities we need. 
One must be very careful, however, since there exist many different parameterizations of the 
experimental phase shifts and of the scattering amplitude so that one must make sure that 
all of the parameterizations are consistent with one another. The second method involves a 
more direct determination of the relevant amplitudes from data. 

a. Forward Amplitudes from Q 2 = 1 — 6 GeV 2 



Wallace has published ]5(| the Pauli amplitudes, A , B°, C°, D°, E° for Hoshizaki's phase 



shift solutions up to Q 2 = 4 GeV 2 . In particular, Wallace uses [54] a more recent solution |58[ 
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than that published in Refs. |55|-|57|) as well as calculating at more energies. Therefore, we 



construct the optical potential strengths in this energy regime by using Eqs. ( |A6[ ) and ( |A7| ) 
for Wallace's amplitudes. 



At higher energies, we use the phase shift analysis of Higuchi and Hoshizaki |59| for pp 
phase shifts at pi a b = 4 GeV, which translates to Q 2 = 5.96 GeV 2 , in order to reconstruct 
the T-matrix and extract the strengths as outlined. Using this method, we were able to 
construct the T-matrix and to reproduce the experimental data and theoretical predictions 
published by Hoshizaki This serves as a consistency check that we are using the phase 
shifts correctly. 

We must stress the fact that at these energies the phase shifts are only loosely constrained, 
at best. Therefore, different phase shift analyses will yield different values for the optical 
potential strengths. However, our strengths are roughly consistent with those found in 
Ref . |5 . 



b. Forward Amplitudes at Q 2 = 9.65 and 20.86 GeV 2 

Phase shift solutions provide a powerful method to summarize the scattering data and to 
predict observables at all momentum transfers. At higher energies, more and more partial 
waves become important (and the data becomes more and more scarce) and phase shift 
analyses become difficult. However, phase shift analyses do exist fl60| , |6"Tf for (lab) beam mo- 



menta of 6 GeV/ c and 12 GeV/ c; in quasi-elastic kinematics, these beam momenta translate 
to Q 2 = 9.65 and 20.86 GeV 2 . However, the lack of constraining data as well as the idea 
that the high angular momentum partial waves are not negligible and must be modelled 
leads us away from the phase shifts. We will, however, take certain qualitative features from 
these analyses; see below. Further, the phase shift analyses do not emphasize the forward 
direction observables over the other angles. Indeed, 

(A9) 
(A10) 

B°, (All) 



A° = 


8tt V 


-ia f ) , 


B° = 


— 167T 


F 2 


E° = 


Aa L 

-167T 


, F 3 
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where atf is the ratio of the real part to the imaginary part of the forward scattering ampli- 
tude and k is the center-of-mass momentum. The other real parts are F 2 = 2 p i Re0 2 (O°) and 
F3 — ^pRe [0i(O°) — 3 (O°)]. The longitudinal cross section, Aa L , is given by the difference 
of measured cross sections with beam and target proton spins aligned and anti-aligned with 
the beam direction with that of the both spins aligned along the beam direction. Similarly, 
the transverse cross section is measured by subtracting cross sections where the beam proton 
spin is aligned with the normal to the beam direction and the target spin is anti-aligned 
minus the cross section where both beam and target spins are aligned with the normal to 
the beam direction. The quantities F2 and F3 are evaluated using dispersion relations [|S3 



Thus, the only required quantity not available directly from the forward cross section 
data is the parameter C°. However, this parameter can be obtained if enough data on the 
spin observables at small —t is known. For instance, in the four component notation of 
Refs. |8^-|86|; we can write two of the spin observables as 

k 2 dn 

-^P n = Rea*e, (A12) 
7r at 

k 2 da 

— A 00s k = —Re a*d sin 9 — Im d*e cos 9, (A13) 

Ti dt 

where 9 is the center-of-mass scattering angle and 

B ^ < A15 > 

c -~h (A16) 
D= - a+b+c - d 

Mk 



E = - a + b + c + d . (A18) 
Aikq 2 



Notice that both of these quantities, P n and A s k (A s k = Aoosfc), vanish in the forward 
direction (9 — e — 0). Their slopes (as a function of q, not —t) at q — 0, however, are 
finite. Taking the derivatives with respect to q = |q| = \f^i and evaluating in the forward 
direction, we obtain 

— — l—r 1 } =Rea*-, (A19) 
71 \dtj f=0 V dq J q=Q q 
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k 2 ( da\ ( dA sk \ 1 « . 

— I — I — = --Rea*d-lmd*-, (A20) 

* \ dt Jt=o\ d q J q =o k 1 
where we have used e/q ~ e'(g) in the forward direction. Note that e/q remains finite in this 
limit. We have also changed variables and written the differential cross section in Lorentz 
invariant form. Now, denoting xr and xi as the real and imaginary parts of the complex 
number x, we can use the above two equations to write 

C °>i k 1 fda\ 

Ak 2 Air andn + aidi \ dt 



d ( dP A i ai ( dA A 

R \dq ) q=0 ^ V dq ) 



(A21) 



C ' a i CR Ana i ( rft )t=o ( dq ) ?=0 ' ^ A22 ' ) 

where a superscript °, indicating the forward direction, is to be understood on the quantities 
Cr, Cj, cir, a I} dn, di. 

Thus, with these equations, and some good spin data, we can extract the optical potential 
strengths directly from data. We note that we have tested these equations for the spin 
parameter C on the 4 GeV/ c data of Hoshizaki and they agree with the value calculated 
strictly by using phase shifts. This gives us confidence that these equations are reliable to 
use at higher energies. 

The experimental values for the a tot , otf, Aa^, Aa^, F 2 and F 3 are summarized in 
Refs. []60| , |6i| . In order to complete the analysis described above, we need information on the 
slopes of P n and A s k at small q. In Refs. pT| , |38| the polarization P n has been measured at 
Piab — 6 GeV jc. Better than that, an empirical fit to the data is given, which can be directly 
differentiated. By averaging these two values from the two references, we obtain the result 

(dj\ 
\ dq 

The uncertainty in this quantity is only about 3%; see Refs. [|B7l , j38|l . At pi a b = 12 GeV/c, the 
polarization has been measured and again an empirical fit is given. Thus, differentiating 
this quantity directly we find that 

fdR 



= 0.494 GeV 1 for p Ub = 6GeV/c. (A23) 

g=0 



dq 



= 0.215 GeV' 1 for p lah = 12 GeV/c. (A24) 

<7=0 
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To proceed, we now notice something amazing. Looking at the data for A s k at pi a b = 
4 GeV/ c we see that the experimental slope of A s j, is zero to a very good approximation. 
Indeed, explicit calculation based on Hoshizaki's phase shifts yield a slope of 0.014 GeV^ 1 . 



There is very little data for A s k at higher energies, yet the phase shift solutions fl60| , |6ll also 
suggest that the slope of A s k is very small. Based on this evidence, we assume that 

at least for energies below pi a b = 12 GeV/ c. 

Using the central values in the published data yields optical potential strengths which 
are larger in magnitude than the lower energy values. However, by adjusting the parameters 
within their experimental error bars, numbers in better agreement (more consistent) with 
the low energy data can be obtained. It is these numbers which are summarized in Table |. 
There is therefore considerable uncertainty in the last two rows of Table Q. In particular, 
the forward scattering data we use to generate the strengths are displayed in Table [TV]. 
The intermediate Pauli amplitudes obtained in the manner described in this appendix are 
displayed in Table [V[ 

Clearly, the most important of these parameters are A , which contains the information 
about the forward scattering amplitude, and the spin-flip parameter C° which is determined 
from the slopes of the forward spin observables. 

APPENDIX B: DWBA AND CT EIKONAL WAVEFUNCTIONS 
1. DWBA Eikonal Wavefunction 



The expression for the DWBA wavefunction, Eq. (|3. 13|) , can be simplified by performing 



the spin algebra. First, we note that in the coordinate system (B,Z,(p), the cross product 
B x p = —Bpcf) {B = |B| and p = |p|). To reduce the wavefunction to its fundamental 
dependence (linear) on the spin matrices, we define the quantities 



U C (B,Z') 

(Bl) 



2ip 

-ipZ'U so (B,Z' 
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4 +) (B,Z) = | [ Z dZ'U so {B,Z'). 

£ J — oo 



Then, we obtain the final DWBA wavefunction as 



where 



*Sii a (R)=JV 



W W(R) 
<)(R) 



X- 



(B2) 



(B3) 



where 



^(+) = e ipZ e u\ 



(+) 



,(+) 



,(+) 



,(+) 



coscug 1 ^ + icr • <fi sinu;^ 

i 



(+) 



p E + M N + V S -V V 



x 



jcr- Z 
+tr ■ B 



(+) / (+) (+) , • \ • (+) / (+) 

COS UJ 2 [0l>1 z u 2 B + % V ) ~ sm U 2 I ^1 B 



SlUO^ (w)^z — UJ^b + W) + COS UJ2 ' ' (^i^B + ^2,Z 



,(+) 



(+) A ,(+) 



(+) 



d_ 

dZ 
d 



1 

B 



(B4) 



(B5) 



(B6) 
(B7) 



and where we have suppressed the argument R in u%, u>2, V s and V v . 

We note that the wavefunction for incoming boundary conditions is obtained simply by 
taking the complex conjugate of the potentials and changing the limits of integration. Thus, 
'I'p.-sr nas the same form as Eqs. ( [BID an d ( P3|) except that uo^ — > oj\ , where i = 1,2 
and 



u[-\B,Z) = -J z ^-[U:(B,Z')+tpZ'U: o (B,Z')} 



dZ' 



(B8) 



(-) 



;b,z) = -| l°°dz'u: (B,z') 



(B9) 



where Z7* represents the complex conjugate of the quantity U. 

It is also instructive to consider the case where the optical potentials are small, in some 
sense. Then, we can perform an order by order expansion of the DWBA wavefunction, in 
terms of the potentials. This is useful as a practical matter because the size of the first-order 



56 



term of the DWBA is a good indicator of how accurate the OBO CT wavefunction will be, 



see Ref. PP|] . Thus, we proceed by expanding in powers of uj\ and neglecting terms of order 
O(iof). We obtain, 



m) = e ipZ V^ (+) - (R) 



(BIO) 



where the subscript m indicates the order of the expansion of the potentials. Decomposing 
the wavefunction in terms of upper and lower components, 



w 



(+) (R) 



X- 



(Bll) 



We present explicit expressions for the order by order expansion of the DWBA wavefunction. 
We only need the first few terms. 



u, 



w, 



(+) 

0,p 

(+) 

0,p 



a p 



E + M N + V S -V V 



(B12) 
(B13) 



„(+) 



w 



(+) 
i.p 



(oui + iu^ a ■ (p 
1 



E + M N + V S -V V 

x\(T-Z 



+<T B 



(B14) 



(B15) 



w 



(+) 

2.p 

(+) 
2,P 



(+), ,(+) 



1 1 r/> , - 2 



£ + Mjv + y s - v v 



a Z 



-f(4 +) ) 



• (+) (+) , • (+) (+) , ■ (+) (+) , • (+) (+) , M ^2 



5 



+ <T B 



(+) (+) ■ (+) (+) , • (+) (+) 

PU0\ UJ 2 1^1 W]b + 2^2 ^2,_B 



,(+), ,(+) ,(+), ,(+) 

"Zu/2 z — 1 2 Z 



(B16) 



(B17) 
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2. CT Eikonal Wavefunctions 



Here, we show the explicit CT eikonal wavefunctions in our various approximation 
schemes. More notation is needed, so we introduce the quantity 

W£?(B, Z, Z>) = U^\B, Z, Z>) - i P Z'u£>\B, Z, Z'), (B18) 

where p/p2j is set to unity. This CT potential (a number) with three arguments is not the 
operator U c ( so y, although there is an intimate relationship. In fact, 

W c ( gJ)(B, Z, Z') = e i{p2m - p)z e - i(pa *- p)z ' 

CO 

x ^(2m|W c(so) (B,Z , )|2j). (B19) 

m=0 

In first-order calculations, a common simple quantity is the integral of these CT potentials 
which connect to the nucleon. Thus, we define the CT integrals as 



^ +) (B, Z) = J ^ ^fflCB, Z, Z% (B20) 
nS +) (B, Z) = | [ Z dZ'u£\B, Z, Z% (B21) 



and their derivatives as 



= (B22) 
^S(R) = 4 fi ! +) ( R )' (B23) 



where i = 1,2. 



a. Order By Order Wavefunction 
The (outgoing) OBO wavefunction is defined such that 

= e ^ E (B24) 

where the index m indicates the order (0-2) of the expansion of the path ordered exponential. 
Using standard techniques and previously defined notation, the result is 
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where 



m OBO . (R) = Af 

m,P,-s_ R V- L V ■ /v 



W 



,OBO(-D\ \ 
OBO 



u, 



OBO 
0,p 



OBO 



cr • p 



°' p E + M N ' 



3 p°(r) y 



X- 



(B25) 



(B26) 
(B27) 



u 



w 



OBO 
1,P 

OBO 
1,P 



(fix + if2 2 cr • 

f <7 Z 



+ 



E + M N 
cr ■ B 



.a 



E + M J 



pQi + i— - iQ ljZ + iQ 2 ,B 
[pfl 2 -i^i B -ifi 2 ,z] >, 



A' 



it 



OBO 
2,P 

OBO 
2,P 



r <r • z 

cr B 



pAi + i 



B 



—%A\ z + iA 



2,B 



+ 



[pAi - Mi,b - iA 2 ,z] \, 



E + M N 

where we have defined two new (outgoing) "second-order functions" by 
r z dZ' r z ' dZ" 



Ai(B,Z) 



x 



(B28) 



(B29) 



(B30) 



(B31) 



—00 2ip J -00 2ip 

z')u^\z") + 5YW£>(Z, z')c/£ ' 0) (^) 

+e ,( P2 - p)( z'-z") ( W (2) (Z) zyjm^jT) + bW*\z, Z')Uff\Z")) 
+e i{ Pi -p)(z>-z") u ^ z ^ z')U^ \Z")\, (B32) 



iA ^ = - B »L2fJ- 



z dZ' r z ' dZ" 
2ip J-00 2ip 



x 



^')f/io°' 0) (z") + wff (z, (z") 

Ap2 - p) (z>-z») (7A2)(z,Z')u£ \Z")+u£\Z,Z')Ui 2 J\Z' 



+e 



:(p 4 - P )(z'-z' 



(B33) 



and we have suppressed the impact parameter inside the arguments of the CT functions for 
clarity. 
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There are a few analytic checks of these expressions. First, in the CT limit, when 
p A ps p 2 ps p becomes very large, then ^1(2) and .4.1(2) approach zero and CT is obtained. 
The next thing to notice is that the OBO\ wavefunction has the same form as the DWBAi 
with the simple substitution u — > Q. Thus, at low energies (below excited state threshhold), 
P2 and P4 are purely imaginary so that U^ so \ — > Uf^l. Thus, insofar as we can take 
Uc(so) = U c (so)> see Eqs. ( |4.11|) and (|4.12| ), we see that f2 — > uo and the DWBAi is obtained, 
modulo, of course, the small CTD effects described above. The same analysis works with 
the second-order calculation. In particular, below threshold, the exponentials damp most 
of the terms so that the only ones which contribute are U^h 80 \ — > U^f) BO y Thus, in this 
low energy limit, A\ — > ujf/2 — uj\ /2 and A2 — > ^1^2- Some simple algebra then yields the 
desired result that, in the low energy limit, the OBO wavefunction approaches the DWBA. 



b. Exponential Approximation Wavefunction 

In this approximation, the CT wavefunction is obtained by taking the first order result 
and exponentiating. Thus, the EA has the same relationship to the OBO\ as does the 
DWBA to the DWBAi. This is equivalent to neglecting the path ordering. Then the 
(outgoing) EA wavefunction can be written down immediately: 



(B34) 



where 



and 



u 



w 



EA, 
P 

EA 



R) 

;r) 



x- 



(B35) 



U EA = e i pZfPl 



W 



EA 



cos 0,2 + icr ■ <j> sin Q2 , 

————{ <j ■ Z cos ^ 2 (fii.z - ^2.b + ip) 
hi + Mn 1 L 



-sinfi 2 [&i,B + ^2,z 



1 

B 

+<x • B sinf2 2 {£li,z - ^2,b + ip) 

+ COsf2 2 (fii,B + ^2,z) 



(B36) 



(B37) 
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We again explore the high- and low-energy limits to perform a check on this wavefunction. 
At very high energies, Huty become very small and approach zero. When this happens, it 
is easy to see that CT is obtained. At very low energies, aside from small quadratic CTD 
effects, the EA approaches the DWBA. 



c. Low Energy Wavefunction 

The CT wavefunction in the LEE approximation is defined so that the expectation value 
of the CT operators in the nucleon is separated out and the rest is treated in perturbation 



theory. In Ref. |3(J we defined the LEE approximation so that the wavefunction reduced 
to the DWBA in the zero'th order. Due to quadratic and CTD effects, this will no longer 
be exactly true. It is necessary to define more notation to obtain the (outgoing) LEE 



wavefunction. Thus we define : 



, n ^ f z dZ' 

J — oo Z,tp 



U^°\B,Z') 
-ipZ>u£°\B,Z') 



(B38) 



u 02 (B, Z) = -j dZ'Uff\B, Z'), (B39) 

L J — oo 

along with the reduced CT integrals 

Q om (B,Z) = n m (B,Z)-cu om (B,Z). (B40) 
With these preliminaries we obtain the (outgoing) LEE wavefunction as 

*S,P,-i*(R) = ^ E *mfp!-sji(^)' (B41) 
m 

where the index m indicates the order of the expansion of the perturbation; m = or 1. 
The result is 

C-JR)^ *^ \X-s R , (B42) 



where 
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„.LEE _ woi 

u o,p - e 
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X 
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(B43) 
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ap E + M N 



cr • Z cos uj 02 [(p - iuj i,z + ^02,b) &i + 
i (vo ljB + ^02,z + j^J B 2 - i&i,z + iB 2 ,b\ 
cr • Z sin uj 02 [{-p + iuj i,z ~ ^02,b) B 2 + 
i (uoi,b + ^02,z + j^J &i + iBi, B + iB 2 ,z] 



cr • Bcoscj 02 [(p - ivoi,z + iu 02 ,B) B 2 - 
i (uoi,b + ^02,z) Bi - iB hB - iB 2 , z ] 
cr ■ B sin uj 02 [(p - iu i,z + ^02,b) Bi + 



i (w i,B + ^02,z) #2 - i#i,z + iB 2 , B ] I , (B48) 
where we have defined two new functions, in analogy to Aip), to be 

^ i(B ' z) = L % L %{ ei{p2 ~ p){z '~ z ' ,) mu -°^ z '^\ z ") 

+ B 2 p 2 el ( P2 - P )(Z>-Z») (0,0) Wi (2) (z> Z ') C/ (2,0) (Z // ) 

+e i(P4- P )(Z'-Z») (0,0) W (4) (Z) Z ')^) (z „ ) | > (B4Q) 

<ft(B, Z) = -Bp g d ^\f^- p){Z '- Zn) m ^l\Z, Z')U^\Z") 
+e i(r2-p)(z>-z») (o,o) w (2) (Z) z')Ug°\Z") 

+e i<p*-p){Z>-Z») (0,0) W (4) (Z) Z ')^) (z „ ) |_ (B5Q) 

We have also defined the new functions 

{m u^l o) {z,z')=u cso{so) {z,z') 

-^- ri(Z - Z,) C(2'). (B51) 

The behavior of this LEE wavefunction at low energies has already been briefly discussed. 
Note that at very low energies, Q 0i — > as well as (°' ^U CSO ( SO ) — > 0. Therefore, except for 
quadratic and CTD effects, the LEE wavefunction approaches the DWBA wavefunction 
for low energies. One can also check this wavefunction at high energies where the potentials 
Q 0i — > — u 0i . Now, this wavefunction claims to be correct only to second order in the CT 
potentials. Thus, if we make the replacement that Q 0i — > — u 0i , we can expand the resulting 
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wavefunction to first order in the potentials cjo«- Then, in addition, letting coscj 2 — ► 
1 — ujq\/2, sinc^02 ~~ * ^02 an d e ^ 01 — ¥ 1 + ^01 + <^oi/2 it is relatively simple to show that, to 
within terms of O(uj 3 ), the LEEi wavefunction approaches the plane wave result for very 
large momenta. In case the reader wishes to verify the above statement, an intermediate 
result is that, in this limit, the B\ = {uiq\ — ^02) /2 + 0(uq) and B 2 = ou 01 u 02 + 0{oJq). 
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FIGURES 

FIG. 1. Schematic drawing of (e,e'p) reaction. 
FIG. 2. Coordinate system used to describe the (e, e'N) reaction. 
FIG. 3. The second-order cross section ratios; M<i = lA4GeV. 

FIG. 4. The second-order cross section ratios; M2 = 1.80 GeV. 

FIG. 5. Total cross section ratios for 12 C. Solid (circles) is DWBA, dotdashed (diamonds) 
includes CT, via EA, for M 2 = 1.44 GeV, dashed (boxes) have M 2 = 1.80 GeV. 

FIG. 6. Total cross section ratios for °Ca. The curves are as in Figure ^. 

FIG. 7. Total cross section ratios for 208 Pb. The curves are as in Figure & 

FIG. 8. Ratio of integrated longitudinal to transverse response for C. Dotted (fancy boxes) 
is plane wave result, solid (circles) is DWBA, dotdashed (diamonds) and dashes (boxes) are 
CT-included cases for light and heavy excited state masses. 

FIG. 9. Ratio of integrated longitudinal to transverse response for 40 Ca. The curves are as in 
Figure |8|. 

FIG. 10. Ratio of integrated longitudinal to transverse response for 208 Pb. The curves are as 
in Figure |8[ 

FIG. 11. Current conservation violations for C. The curves are as in Figure ||[ 
FIG. 12. Current conservation violations for 40 Ca. The curves are as in Figure ||. 
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FIG. 13. Current conservation violations for 20S Pb. The curves are as in Figure ||. 

FIG. 14. Differential unpolarized cross section and normal polarization for 12 C at 
Q 2 = 0.96 GeV 2 . The curves are as in Figure ^. 

FIG. 15. Differential unpolarized cross section and normal polarization for 12 C at 
Q 2 = 20.86 GeV 2 . The curves are as in Figure |8|. 

FIG. 16. Effect of different electron kinematics at Q 2 = 5.96 GeV 2 . The dotdashed curve is the 
EA for M2 = 1.44 GeV, the dashed curve is the OBO\, the dotted curve is the LEE\. The solid 
curves are the DWBA and the plane wave result, the context identifies them. 

FIG. 17. Differential unpolarized response functions for 208 Pb at Q 2 = 3.25 GeV 2 . The curves 
are as in Figure ||. 

FIG. 18. Differential normal response functions for 208 Pb at Q 2 = 3.25 GeV 2 . The curves are 
as in Figure ||. 

FIG. 19. Cross section ratios for finite-sized wavepacket. 
FIG. 20. Current conservation violations for finite-sized wavepacket. 

FIG. 21. Effect of Fermi motion. 

FIG. 22. Total cross section ratio for 12 C with no lower components. The upper CT curve and 
the lower DWBA curve are calculated with wavefunctions having no lower components. The point 
is that the ratio is independent of the lower components. 

FIG. 23. Current conservation violations for 12 C with no lower components. In this figure, 
the dotted curve is the Born calculation with no lower components; the other curves are as in 
Figure [2^. 
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FIG. 24. Ratio of integrated longitudinal to transverse response for 12 C with no lower compo- 
nents. The curves are as in Figure W\. 
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TABLES 



TABLE I. Energy dependent strengths of optical potentials. Q 2 is in GeV 2 and pi a b is in GeV. 
r = -4iripi ab /M N . 



Q 2 


Plab 


rF°(fm 2 ) 


rF°(fm 2 ) 


Po rF°(MeV) 


Po rF°(MeV) 


0.96 


1.1 


-9.566+ 1.768 * 


5.720-2.291*' 


-313.3+ 57.91*' 


187.4 - 75.05 i 


1.88 


1.7 


-9.550+ 3.748 * 


4.970-3.895*' 


-312.8+122.8* 


162.8 -127.6*' 


2.38 


2.0 


-10.46 + 4.219* 


4.912-3.936* 


-342.6+138.2 i 


160.9 -128.9*' 


3.25 


2.5 


-7.243+ 4.047 * 


3.396-3.573*' 


-237.3+132.6*' 


111.2 -117.0*' 


4.14 


3.0 


-10.75 + 6.735 * 


3.836-4.124* 


-352.0+220.6*' 


125.6 -135.1*' 


5.96 


4.0 


-7.091+ 2.147* 


2.155-2.540*' 


-232.3+ 70.33*' 


70.59- 83.20 i 


9.65 


6.0 


-7.681+ 2.398*' 


1.831-2.380* 


-251.6+ 78.56*' 


59.98- 77.96 i 


20.86 


12.0 


-7.631+10.73 i 


1.168-2.809* 


-250.0+351.4* 


38.27- 92.03 i 


TABLE II. Cross section and ratio of real to imaginary parts of forward scattering amplitude, 
as given by optical potential strengths. 


Q 2 

(GeV 2 ) 




Elab 

(GeV) 




a 
(mb) 


Of 

(unitless) 


0.96 




1.45 




30.1 


0.43 


1.88 




1.94 




47.7 


-0.17 


2.38 




2.21 




47.5 


-0.22 


3.25 




2.67 




46.0 


-0.40 


4.14 




3.14 




44.3 


-0.30 


5.96 




4.11 




42.2 


-0.26 


9.65 




6.07 




40.8 


-0.32 


20.86 




12.04 




39.6 


-0.29 
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TABLE III. Calculational kinematics. 



Q 2 

(GeV 2 ) 


Ei 
(GeV) 


0e 

(deg) 


( da ) 

\ dn k' J Mott 
(nb) 


V L 


V T 


Vtt 


V LT 


0.96 


4.0 


15.1 


1091. 


0.62 


0.41 


0.39 


-0.70 


1.88 


4.0 


22.8 


203.0 


0.43 


0.37 


0.33 


-0.54 


2.38 


4.0 


27.0 


103.6 


0.36 


0.36 


0.30 


-0.48 


3.25 


4.0 


34.8 


36.90 


0.27 


0.36 


0.26 


-0.41 


4.14 


4.0 


44.6 


13.40 


0.21 


0.40 


0.23 


-0.36 


5.96 


4.0 


84.5 


0.8865 


0.14 


1.01 


0.19 


-0.41 


5.96 


6.0 


34.5 


17.07 


0.14 


0.28 


0.19 


-0.25 


5.96 


11.0 


22.7 


140.7 


0.14 


0.20 


0.19 


-0.24 


9.65 


15.0 


14.7 


85.35 


0.072 


0.15 


0.13 


-0.14 


20.86 


21.0 


18.2 


18.24 


0.021 


0.098 


0.072 


-0.060 



TABLE IV. Forward scattering data at pi a b = 6 and 12 GeV/ c. 

Piab = QGeV/c 

a tot = 40.75 mb a f = -0.32 

Aa T = 0.35 mb F 2 = -4.60 GeV' 1 

Aa L = -1.04 mb F 3 = 4.60 GeV' 1 
(wLo = 93.0 mb (^) 



Piab = 12 GeV/ c 



0.494 GeV- 



i 



a tot = 39.60 mb a f = -0.29 

Act t = 0.01 mb F 2 = -6.05 GeV' 1 

Aa L = -0.73 mb F 3 = 3.55 GeV' 1 

(w) t=0 " 65.0 mb - 0.215 Gey 



-i 
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TABLE V. Pauli amplitudes for pp scattering 



Q 2 


A 


B° 


C° 


£° 


(GeV 2 ) 


(GeV- 2 ) 


(GeV- 2 ) 






0.96 


3.08 - 1.31 z 


-0.22 + 1.21 * 


-1.76 - 6.55 z 


0.69 + 0.02 * 


1.88 


4 87 + o 81 * 

A • KJ 1 | KJ * KJ -L V 


-0 23 + 1 03 * 

\J - ".7 | J. • \J KJ v 


-3 13 - 5 44 % 

KJ * a kj kj • a jr i 


82 - 05 * 


2.38 


4 85 + 1 06 * 

A • KJ KJ | -1 • \J KJ V 


-0 38 + 71 * 

KJ • * ' W | KJ - 1 -1. r 


-3.11 - 5.22 i 


81 + 16* 

\J * KJ -L. \ KJ • A \J L/ 


3.25 


4.70 + 1.86 * 


-0.14 + 0.48*' 


-2.81 — 3.35 i 


0.29 + 0.16*' 


4.14 


4.53 + 1.34*' 


-0.05 + 0.47*' 


-3.41 - 3.95 i 


0.15 - 0.01* 


5.96 


4.31 + 1.13* 


-0.04 + 0.40*' 


-1.62 - 2.13? 


0.12 - 0.09*' 


9.65 


4.16 + 1.33* 


-0.02 + 0.19*' 


-1.34 - 1.62* 


0.07- 0.00*' 


20.86 


4.05 + 1.17* 


-0.00 + 0.13*' 


-1.61 - 0.87*' 


0.04- 0.05*' 
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